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EXECUTIVE SUMMARY 


EXECUTIVE SUMMARY 


Consistent with the statement of work in the original proposal to NASA, several projects 
were begun in 1988. This report covers the progress on these projects and other 
accomplishments for the period from August 1988 through February 1989. The work 
planned for the period from March 1989 to March 1990 is also described. 

All of the tasks outlined in the proposal were considered during one meeting of the 
NASA Technical Representative Committee and two meetings of the SERC Advisory 
Committee held at the Center. Based on advice from these committees and extensive 
discussions between the Center Director and the two Principal Investigators, those listed 
in the Table of Contents were selected, and their selection formalized through an 
internal summary document listing specific goals, measurable results, and faculty as 
well as student level of effort. 

Though categorized as Propellant Processing, Material Properties, Systems Optimization, 
or Data Base Development, these projects actually represent beginning work in a 
number of distinct but closely related disciplines: (1) propellant and propulsion 

optimization, (2) automation of propellant processing with quantitative simulation, (3) ore 
reduction through chlorination and free radical production, (4) characterization of lunar 
ilmenite and its simulants, (5) carbothermal reduction of ilmenite with special reference 
to microgravity chemical reactor design, (6) gaseous carbonyl extraction and 
purification of ferrous metals, (7) overall energy management, and (8) information 
management for space materials processing. All yielded sufficiently promising results 
to recommend continuation during 1989-90. 

Additionally, low levels of effort were devoted to the development of anhydrous 
cementing and construction materials from lunar soils, and to molten metal-molten salt 
separation research. The first of these has been redirected to feature melt binding 
methods and high purity iron characterization, the second to feature solvent extraction 
with recyclable organic chelating agents. 

The technical highlights of the 1988-89 work are: 

1. Extensive performance calculations for unconventional propellant 
combinations. 

2. Evolution of a quantitative methodology for Figure-of-Merit (FoM) 

calculations for various space mission scenarios. 

3. Design and construction of the process control hardware in one 

propellant processor. 

4. Development of the fundamental criteria for automatic control of 

propellant processors. 

5. Demonstration of significant gains through plasma processing of ores. 

6. Demonstration of the relative merits of carbothermal reduction of ores. 

7. System development for gaseous carbonyl extraction and purification 
research. 

8. Proof-of-principle delineation of pinch technology for energy 

management. 

Details of this work and other developments are provided in the following sections. 
Though still in its initial stage, the program has already made significant progress 
toward the overall goal of exploiting extraterrestrial resources for space missions (see 
Figures 1.1 and 1.2). 


IV 


Several combinations of unconventional propellants were shown to have acceptable 
performance in the light of launch mass reductions. One specific propellant processor 
was acquired and automation procedures started; both hardware and software are under 
development. The systems analysis mission optimization program, even in its present 
simplistic configuration featuring a popular spreadsheet interfaced with an expert 
system, yielded valuable insight into the tradeoffs involved in performance and earth- 
transported resources. The chlorination processing of ores has shown promise in terms 
of rates and compactness of hardware, while pointing up the potential of free radical 
production. The carbothermal task group has not only evolved a flowsheet process 
model, but also built much of the equipment necessary to implement it on a laboratory 
scale. The carbonyl task group has made similar progress and placed orders for the 
required equipment. The "pinch technology" concept evolved by the energy 

management group promises to maximize the use of available energy. 

To support the above and future work of the Center the following specialized items of 
equipment were acquired: (1) a variable ionization mass spectrometer, (2) a Baker- 

Perkins propellant processor, (3) an Inframetrics color-coded video with dedicated 
processor, (4) an IBM-compatible data acquisition and processing system, together with 
appropriate NASA (LeRC), Air Force, and industry propellant performance software. 

A computing and communications network featuring high performance workstations 
(currently SUN 386i, IBM PS2/80, MACII) accessible to all investigators through ethernet 
links was also established at Center Headquarters. Transparency to SUN, VAX, MAX, 
and IBM systems is being stressed. Plans were also completed to install a microwave 
link between the University’s local and national network and the Center’s off-campus 
facility to facilitate more general access. 

Another highlight of the startup program was the large number of distinguished visitors, 
most of whom came and presented seminars at no expense to the Center. Widespread 
faculty interest was generated and several high-quality students joined the program, 
due in no small measure to these seminars. 

Twelve graduate students, including three women and four ethnic minority students, 
were supported during this period, and many applications for future participation were 
received— partly as a result of the seminars, and partly from the wide distribution of a 
flyer and poster describing the Center’s mission. 

Following the recommendation of the NASA Technical Representative Committee, an 
attempt was made to promote cooperation with other NASA SERCs. In terms of 
programmatic emphasis, those Centers that have the closest relation to our own are the 
ones at the Pennsylvania State University (space propulsion), University of Cincinnati 
(monitoring of rocket engine performance), and University of Colorado (space 
construction). A Program Coordination Task Group was initiated and a first meeting 
held at the Center which was attended by representatives of these Centers, as well as 
representatives from the Los Alamos National Laboratory, the Space Sciences Institute, 
and the California Space Institute. 

It is perhaps worth noting that during this initial period of operation, well over two- 
thirds of the total direct-cost funding received from NASA was expended on support of 
the projects listed in the Table of Contents, together with nearly the same proportion of 
the matching funds received from the University of Arizona. Since all faculty 
investigator salary allowances and graduate fellowship stipends are included, 
preserving this ratio is believed to be the most effective way to maximize return on 
investment, both in terms of space processing research and space infrastructure 
development. 
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Several new projects are planned for 1989-90; these include: 


1. Determination of Lunar Ilmenite Abundances From Remotely Sensed 
Data (R. B. Singer, Lunar and Planetary Laboratory) 

2. Continuous Monitoring of the Lunar or Martian Subsurface Using On- 
Board Pattern Recognition and Neural Network Processing of Rover 
Geophysical Data (C. Glass and B. Sternberg, Mining and Geological 
Engineering Department) 

3. Detection of Water on Near-Earth Asteroids (L. A. Lebofsky, Lunar and 
Planetary Laboratory) 

4. Analysis of Near-Earth Asteroids and Cometary Cores (T. Gehrels, 
Lunar and Planetary Laboratory) 

5. Production of Solar Reflectors From Lunar Materials for Solar Power in 
Space (R. J. Thibeault, Science Applications International, San Diego) 


VI 



Figure 1.1. One concrete example. 
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Figure 1.2. Engineering tasks. 



I. PROPELLANT PROCESSING 



Proof-of-Concept Automation of Propellant Processing 
K. Ramohalli and P. Schallhorn 
Department of Aerospace and Mechanical Engineering 
The University of Arizona 
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Abstract 




For space-based propellant production, automation of the process is needed. 
Currently, all phases of terrestrial production have some form of human interaction. A 
mixer has been acquired to help perform the tasks of automation. We have designed, 
built, and installed a heating system to be used with the mixer. Tests performed on the 
heating system verify design criteria. An IBM PS/2 personal computer has been 
acquired for future automation work. It is hoped that by the end of the next academic 
year, the mixing process itself will be automated. This is a concept demonstration 
task-proving that propellant production CAN be automated reliably. 


Introduction 

The research work deals with the autonomous production of propellants. Because 
80% to 90% of a spacecraft’s weight is propellant, it is advantageous to produce 
propellants in strategic locations en route to, and at, the desired mission destination. 
This will then reduce the weight of the spacecraft and the cost of each mission. Since 
one of the primary goals of the space program is safety, a totally automated propellant 
production system is desirable. This system would thereby eliminate the need for the 
constant human intervention currently required in production of many propellants. This 
enables the exploration of space to be more than the search for, and production of, 
propellants. As a proof-of-concept demonstration, one specific case was chosen for 
this study— composite propellant production; the principle is more important than the 
application. 


Background 

Currently, composite solid propellant production is done with constant human 
intervention. Using a control room, man has total control over all aspects of the 
propellant production. This is fine on Earth, but it is too costly in space. Thus, the 
need for automated composite propellant production exists. 


Approach 

We are currently completing testing of a heating system, which was designed by 
the student (Paul Schallhorn), for the one-pint mixer that is to be used for this project. 
Because composite propellant production requires mixing the ingredients at two 
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constant temperatures (160 and 140°F), a self-contained water-heating system is 
required for space-based operation. Such a system is shown in Figure 1.1. This 
system provides the required temperatures and only needs an electric power source to 
drive the pump motor and heat the water heaters. This is not unrealistic considering 
that electricity is also required for the mixer and controlling computer. 



Figure 1.1. The heating system. 


One approach, therefore, is to use a personal computer to control the introduction 
and mixing of the composite propellant ingredients to the mixer (making sure that 
temperature is constant on the walls of the bowl, detecting local "hot spots" within the 
mixture, and taking in-situ measurements of the viscosity of the mixture to check if it is 
within an acceptable range). Then, pump the mixture, via computer programs, into a 
cast which will be placed in an oven for curing and then stored for future use. 

Results to Date 

The major results to date are as follows: 

1. A used Baker-Perkins PX-2 mixer was acquired; this introduces a factor of 6 
cost reduction (see Fig. 1.2 for the complete mixer setup). A heating system 
was required for its operation. 

2. In September 1988, Schallhorn designed the heating system to be used for the 
mixer (see Fig. 1.1). It was determined that the minimum volumetric flow rate 
for the heating system for a 1 -degree temperature drop acrosss the mixer 
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operating at steady state was 2.5 gallons per minute. Therefore, we selected a 
pump with a volumetric flow rate of 4.4 gallons per minute to ensure a 
negligible temperature drop across the mixer bowl. Since only two 

temperatures are needed, it was logical to have two separate reservoirs, each 
at one of the required temperatures. We chose to have both reservoirs be hot 
water heaters. Because we only had 120-volt a/c power available, we had to 
choose the most efficient heater size on the market. As we began to search 
for heaters for the project, it was discovered that the same heating element 
was commonly used in different-sized 120-volt water heaters. This made it 
clear that for maximum water heating, the smaller the water heater, the more 
advantageous. That was the basis for the selection of two 10-gallon water 
heaters (see Fig. 1.3). The system uses distilled water to eliminate the 

possibility of scale buildup in the system. To further ensure the cleanliness of 
the water in the system, a filter is placed in the system immediately following 
the pump (see Fig. 1.4). 

3. Acquisition of the components of the heating system was begun in October 

1988. By the middle of November, all of the components were in and the 
heating system was assembled. 

4. Initial verification of the temperature profile of the heating system was begun 
in December 1988. Verification of the heating system continued through March 

1989, including verification of flow rate and the time required to heat the 
system from a cold start. 

5. In August 1988, research was begun to determine which personal computer to 

purchase for this project. By the end of September, an IBM PS/2, Model 80 
was selected, with an Intel 80386 microprocessor operating at 20 MHz, a 
115-megabyte hard disk drive, and 2 megabytes of RAM. The computer was 
ordered at the end of September, along with the following peripherals: a 

14-inch monitor, a 80387 math coprocessor, a modem, a 5.25-inch external 
diskette drive, additional memory, a mouse, and a Hewlett Packard Laserjet II 
printer. Due to shipping problems from IBM, the computer did not arrive until 
late in January, and the peripherals did not arrive until early February. By the 
middle of February, the computer system was operational. This computer 
system will be used on various other NASA Center projects, also. 

Summary and Future Work 


In summary, this task strongly suggests that there is a need for automated 
production of propellants for space-based propellant production. We have also seen 
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that there is no current system to produce composite propellants without human 
intervention. A mixer has been acquired to help perform this task. We have designed 
and built a heating system to be used in conjunction with the mixer to maintain 
constant mixing temperature. The heating system has been, and is continuing to be, 
tested under operational conditions for design verification. We have acquired an IBM 
PS/2 personal computer for the computer portion of the automation. 

For the 1989-90 academic year, the student plans to begin his Ph.D. research, 
which will consist of the actual automated propellant production. During the year, we 
will begin to automate the mixing process itself. It is hoped to have the computer 
control the addition of each ingredient from a "hopper" (yet to be built) to the mixer at 
required times and have the computer control the mixing of the ingredients for the 
required amount of time. We also plan on building and installing the in-situ viscosity 
measuring device for future integration into the automation system. 

Acknowledgments 

We would like to take this opportunity to thank University of Arizona technician 
Gary Hopkins for his help in assembling the heating system. We would also like to 
acknowledge Richard Wilson, who is doing the automatic controls portion of the 
heating system. 



Extraction of Volatiles and Metals From Extraterrestrial Ores 

J. S. Lewis 

Lunar and Planetary Laboratory 

The University of Arizona j 1 v 

During the start-up phase of this project, we have concentrated upon planning our 
laboratory efforts for the second fiscal year. The principal motivating purpose of our 
work has been to identify extraterrestrial materials, processes, and products associated 
with the production of propellants in space, including the most complete possible 
conversion of the feedstocks for propellant production into useful products with the 
minimum feasible expenditure of energy. We have tentatively identified and are 
beginning laboratory research on several processes that promise very large increases in 
the mass of useful products at the cost of only modest increases in energy 
consumption. Processes for manufacturing propellants then become processes for 
making propellants plus metals and refractories. It is the overall yield of useful 
materials per unit expended energy that matters, not simply the yield of propellants. 
Three tasks have been undertaken to date: 

1 . Literature search and compilation of a dBase 111 data base on space materials 
processing. This work is being done mainly by an undergraduate, Leo 
Masursky, with many of the references being supplied by graduate students 
working for various SERC projects. The data base is roughly half finished. 
Completion of the data base is expected this summer. Thereafter, a low 
level of effort should suffice to keep it current. The data base is designed 
with numerous subject flags to permit efficient searches on any of a large 
number of subjects. This is the heart of a computerized annotated data base 
to be prepared by Dr. Andy Cutler. Andy has assisted us in the design of 
the data base. 

2. Gaseous carbonyl extraction and purification of ferrous metals. This work is 
being carried out principally by a graduate student, Mike Nolan, and 
Professor Lewis. Our activities during the first fiscal year have been 
concentrated on planning and equipment design of our carbonyl apparatus. 
Orders have gone out for a number of crucial components of the gas system. 
Renovations of the lab space have been completed, and the lab is now 
occupied. The gas line is being designed to have the capability of mixing 
several gases in precisely measured proportions, including both carbonyls 
and inert carrier gas (nitrogen or argon); compressing the mixtures to a 
presssure of 10 atmospheres; and regulating temperatures over the range 0 to 
250° C. We have also drawn up plans for extraction experiments on native 
meteoritic metal alloys, magnetic separates from lunar fines, cathode metal 
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deposits from the lunar magma electrolysis scheme (Haskin), and the 
byproduct elemental iron from reduction of lunar ilmenite (Shadman). 

3. Characterization of lunar ilmenite and its simulants. This work is the principal 
responsibility of a graduate student, Melinda Hutson. During the past year, 
she has completed a literature search on the chemical and physical 
properties of lunar and meteoritic ilmenite and of schemes for beneficiation 
of ilmenite from lunar regolith (see Appendix A). Certain achondritic 
(basaltic) meteorites contain either large crystals or large concentrations (not 
both) of ilmenite which has a chemical composition similar to, but not 
identical with, lunar ilmenite. The meteoritic ilmenite is free of ferric iron, 
the component of terrestrial ilmenite whose presence invalidates its use as a 
lunar simulant. Museum catalogs have been examined to identify the best 
sources for the most promising meteorites, and samples are on order. 

Experiments will be conducted by Shadman on lunar, synthetic, and 

meteoritic ilmenite to determine which, if any, are satisfactory kinetic and 
thermodynamic simulants of real lunar ilmenite. With Joaquin Ruiz, we are 
preparing a proposal for Apollo 17 lunar regolith samples to be used in 

studying the physical properties of lunar ilmenite and in designing and 

testing beneficiation schemes for extraction of ilmenite from typical lunar 
feedstocks. 
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Production of Oxygen From Lunar Ilmenite 
F. Shadman and Y. Zhao 
Department of Chemical Engineering 
The University of Arizona 



Abstract 

The overall objective of this project is to develop a novel carbothermal reduction 
process for production of oxygen from lunar ilmenite. The specific objective is to use 
a reaction sequence in which a wide variety of carbonaceous compounds (including 
carbonaceous wastes) can be used as reducing agents. During the first phase of this 
work, two reactor systems were designed, constructed, and operated to investigate the 
fundamentals of the reactions important in this process. One system is a small fluidized 
bed, and the other is a thermo-gravimetric reactor system. Preliminary experiments on 
synthetic ilmenite are conducted to study the effect of carbon type, carbon loading, 
temperature, and gas flow rate. Results indicate that a reaction path based on carbon 
gasification can be employed to promote the overall kinetics. A unique temperature 
and concentration-programmed reaction procedure is being developed for rapid 
parametric study of the process. 


Introduction 

Oxygen is a consumable material which needs to be produced continuously in most 
space missions. Its use for propulsion, as well as life support, makes oxygen one of 
the largest volume chemicals to be produced in space. Production of oxygen from 
lunar materials is of particular interest and is a very attractive possibility. 

There is a significant amount of oxygen on the moon, although very little is readily 
available as water or gaseous oxygen. The only practical source of oxygen on the 
lunar surface is igneous materials, which contain typically 40-50% oxygen as oxides. 
Minerals present in these rocks include ilmenite (the most abundant opaque material in 
lunar rocks), anorthite, and olivine. The igneous rocks have been pre-crushed to form 
regolith, the lunar equivalent of soil, by meteoritic bombardments. This reduces the 
need for crushing and simplifies the mining and separation problem greatly . 1 ’ 2 

Extraction of oxygen from iron oxide present in lunar ilmenite is of particular 
interest because it is energetically more favorable than extraction from silicon, 
aluminum, titanium, calcium, or magnesium oxides. Iron oxide reduction is also 
attractive because of its potential for producting iron as a co-product. A number of 
processes have been suggested for oxygen production from ilmenite in lunar regolith. 
Table 1.1 shows a list of these processes. Most of these processes require imported 



1-11 


reagents to be recycled in the process. Therefore, low loss per mass is a critical 
factor and requirement for process feasibility. 


Table 1.1. Selected proposed processes. 


Hydrogen reduction 

Low-temperature carbothermal reduction 
Slagging (high-temperature) carbothermal reduction 
Slag electrolysis 
Leaching: Fluoroacid 
HCI 

Caustic 
Carbo-ch I or i nat i on 

Vaporization/rapid quenching, arc electrolysis, plasma reduction, alkali metal thermal 
cycling, vapor-phase pyrolysis, etc. 


Among the proposed processes, the hydrogen and carbothermal reductions of 
ilmenite appear very promising. Hydrogen reduction has a relatively simple process 
configuration; the individual steps are relatively well studied. However, the major 
problem is the large heating and cooling loads required to condense the water and 
then heat the hydrogen to its reaction temperature. Handling and storage of large 
amounts of hydrogen is also a problem. Hydrogen reduction of lunar ilmenite has been 
studied by several investigators. 3 Carbothermal reduction of terrestrial ilmenite has also 
been extensively studied, but the application to lunar ilmenite is still an open area for 
investigation. 4 

Even though a considerable body of literature has been published on various 
aspects of the reduction of terrestrial ilmenite, the mechanism of transformation and the 
rate-controlling step in the overall reduction have not been clearly established. Most 
studies on terrestrial ilmenite show that both hydrogen reduction reactions take place in 
two main stages. In the first stage, the ferric iron component of the ilmenite was 
rapidly reduced to ferrous iron. In the second stage, which was much slower than the 
first stage, the ferrous iron was reduced to the metallic state. During the carbothermal 
process, terrestrial ilmenite starts to react at about 860° C at the contact points between 
the reactants. Up to 1020°C, solid-state reduction appears to be the main reaction 
mechanism, while above this temperature, a rate increase has been observed and has 
been contributed to a change in the mechanism to a gaseous reduction of ilmenite by 
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regenerated CO. 5 - 6 It has been reported that the reaction rate between synthetic 
ilmenite and graphite is increased significantly by the addition of ferric chloride, which 
promotes the nucleation of iron, and is decreased by the addition of rutile. 

The above mechanisms and the available kinetic information are not directly 
applicable to lunar ilmenite, which is very different from terrestrial ilmenite in both 
physical and chemical properties. Moreover, the emphasis in reduction of terrestrial 
ilmenite has been on metal oxide recovery and not on oxygen production. However, 
the existing technical background on carbothermal processing of terrestrial ilmenite 
indicates that its application to lunar ilmenite should be pursued. In comparison with 
the hydrogen reduction, the carbothermal process has a more complicated flow 
configuration, but is energetically more favorable, and can potentially use a wide 
variety of carbonaceous compounds as reducing agents. A comparison of the primary 
features of the carbothermal and hydrogen reduction processes can be made using 
Tables 1.2 and 1.3. 


Table 1.2. Hydrogen reduction. 

Advantages 

Simple flow sheet 

Relatively well-known chemistry 

Continuous operation 

Disadvantages 

Low hydrogen utilization 

Problems in hydrogen transport, handling, and storage 

High heating/cooling demand 

No useful co-products (low integration value) 

Key Research Areas 

High-temperature water dissociation 
High-temperature ceramic membranes 

Sensitivity to the physical and chemical properties of the feed 
Process design 
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Table 1.3. Carbothermal reduction. 

Advantages 

Known technology for terrestrial ilmenite 
Available carbonaceous makeup 
Almost 100% utilization of carbon 
Useful co-products (high integration value) 

Ease of carbon transport, storage, and handling 

Disadvantages 

More complex flow sheet 

More solid handling in continuous operation 

Key Research Areas 

Reaction kinetics and mechanism for lunar ilmenite 
Ceramic membranes application 
C0 2 reverse fuel ceil 

Sensitivity to the physical and chemical properties of the feed 
Process Design 


Objectives 

The overall objectives of this study can be described as follows: 

• Study of the mechanism and kinetics of carbothermal reduction of simulated 
lunar ilmenite; determination of the rate-limiting steps; investigation of the 
effect of impurities, particularly magnesium; search for catalysts suitable for 
enhancement of the rate-limiting steps. 

• Development of new reduction methods which are based on the use of a waste 
carbonaceous compound as a carbon source for the process; development of a 
novel carbothermal reaction path which utilizes gasification of a carbonaceous 
solid for reducing gaseous species (hydrocarbons and carbon monoxide) to 
facilitate the reduction reaction kinetics and make the process more flexible in 
using various forms of carbonaceous feeds. 

• Development of advanced gas separation techniques, including the use of high- 
temperature ceramic membranes. 

• Development of an optimized process flow sheet for carbothermal reduction 
and comparison of this process with the hydrogen reduction scheme, as well 
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as a general comparison with other leading oxygen production schemes; use of 
new and advanced material processing and separation techniques. 

Method of Approach 

In the first phase of the work, the emphasis has been on the development of 
experimental methods for studying the kinetics and mechanism of reduction of 
simulated lunar ilmenite. The experimental facilities and procedures are summarized 
below. 

Experimental Apparatus. Two reactor systems have been designed, fabricated, and 
put into operation during the past 9 months: A mini-fluidized bed reactor for quick 

response times and a thermo-gravimetric reactor system with an electronic 
microbalance for the continuous measurement of rate and sample conversion. 
Schematic diagrams of these two reactor systems are given in Figures 1.5, 1.6, and 1.7. 
The components of these two reactor systems are described below. 



Figure 1.5. Mini-fluidized bed reactor system. 
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Figure 1.6. Micro-gravimetric reactor system. 



Figure 1.7. Schematic of the quartz reactor. 
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a. Gas Preparation Section: This includes a UHP nitrogen delivery system with 
critical flow elements for mass flow control. 

b. The Mini-Fluidized Bed Reactor: This is a quartz reactor with a U 

configuration. The reactor is mounted on the stand and can be moved 
vertically in and out of an electric furnace. A chromel/alumel type K 
thermocouple is placed adjacent to the fluidized bed to monitor the 
temperatures independently during the course of an experiment. The ilmenite 
particles are supported by a stainless steel screen fitted inside the reactor. 
This reactor has low dead volume and a very rapid response. 

c. Thermo-Gravimetric Reactor: In this reactor system, the sample is suspended 
on one arm of an electronic recording balance. A specially designed quartz 
reactor, which surrounds the sample with ports for inlet and dilution gas 
entry, is mounted vertically on wheels. This configuration allowed rapid 
start-up and quenching of the experiments for accurate measurement of initial 
and instantaneous rates. The furnace is equipped with a PID controller for 
both steady and temperature-programmed reactions. 

d. Gas Analysis Section: This includes nondispersive infrared analyzers for 

measurement of carbon monoxide and carbon dioxide, as well as a gas 
chromatograph for a more complete gas analysis. 

e. Other facilities used in our experimental study include an X-ray diffraction 
spectrometer, Mossbauer spectrometer, and scanning electron microscope 
with EDX. 

Experimental Procedure. The work so far has focused on the study of the 
mechanism and kinetics of reactions important for carbothermal reduction of lunar 
ilmenite. The samples used have been synthetic ilmenite (iron titanate) obtained from 
CERAC, Inc., with a purity of 99.9% and particle size between 63 and 74 microns. To 
simulate lunar ilmenite, it was important to make sure that the iron in the sample was 
only Fe ++ . To confirm this, Mossbauer spectroscopy was performed on the synthetic 
ilmenite. A typical Mossbauer spectra of the sample is shown in Figure 1.8. The 
results indicate that the amount of impurities, including Fe +++ , was less than 0.8% on a 
metal basis. 

In each experimental run, a well-characterized mixture of ilmenite and carbon with 
known particle size and composition was placed in the reactor while the furnace was 
heated to a desired temperature. The gas flow inside the reactor was then stabilized 
and maintained for purging the reactor and the sample while the furnace temperature 
approached the pre-set value. To start the run, the reactor was lowered into the furnace 
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Figure 1.8. Mossbauer velocity spectrum for synthetic ilmenite particles. 
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(in the fluidized bed system) or the furnace was raised around the reactor (in the 
thermo-gravimetric system). The course of the reaction was studied by monitoring the 
gas composition and the sample weight in the gravimetric reactor. 

The other methods used for characterization of the starting ilmenite were scanning 
electron microscopy and X-ray diffraction spectroscopy. Application of these 
techniques to samples at different levels of conversion can help us to understand the 
morphological and phase changes that occur in the samples during the course of 
reaction. 

Discussion of Results 

Due to the heterogeneous nature of the reduction reactions, the kinetics of the 
process depend on the physical structure of the solid. Figures 1.9 and 1.10 are the 
micrographs of the original synthetic ilmenite. The small particles on the surface are 
tiny chips generated during the grinding process. Figures 1.11 and 1.12 show partially 
reacted samples and, in particular, the porosity developed during the reaction. The 
SEM micrographs show a significant change in the structure and porosity of the 
particles as the reduction reaction proceeds. It appears that the new products formed 
on the surface of the particles have a lower melting point compared with that of 
ilmenite (melting point of iron titanate is 1327°C). Further characterization of this phase 
change is in progress. 

So far, two types of experiments have been performed in the reactor systems. 
Most of the results are tests of the accuracy, reproducibility, and calibration of the 
newly developed reactor systems. More recently, however, experiments for determining 
the kinetics and mechanisms have been started. Among the first type of experiments 
are tests of the differential nature of the reactors. This is important to make sure that 
the kinetics are characteristics of the reaction and not the reactor geometry; this is an 
important point often neglected in similar studies. As an example, tests with different 
sample mass were conducted and the initial rates were compared. The results are 
shown in Table 1.4. The fact that the initial rate did not change significantly with 
sample size confirms that the reactor is differential and that the rates are intrinsic 
values. Similar results were obtained when flow rate was varied. 

One of the important process parameters is the optimum carbon/ilmenite ratio in 
the reactor. The preliminary results are shown in Figure 1.13. It appears that the initial 
rate increases slightly with an increase in the carbon content. The significant 
difference is in the rate profile. As carbon content decreases, the variations in rate 
with time and conversion become more significant. The variations in rate are very 
important in the design of a large-scale reactor. 
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Figure 1.9. The surface microstructure of synthetic ilmenite (particle size ss 70 /zm; 

sample coated with Au/pd); specimen tilt, 30°; magnification, 600X; left bar 
length, 10 /zm. 



Figure 1.10. The surface microstructure of synthetic ilmenite (particle size ^ 70 /zm; 

sample coated with Au/pd); specimen tilt, 30°; magnification, 1030X; left 
bar length, 10 /zm. 
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Figure 1.11. The surface microstructure of the reduced synthetic ilmenite (particle size 
70 /im; sample coated with Au/pd); specimen tilt, 30°; magnification, 
490X; left bar length, 10 /im. 



Figure 1.12. The surface microstructure of the reduced synthetic ilmenite (particle size 
2 70 /an; sample coated with Au/pd); specimen tilt, 30°; magnification, 
1050X; left bar length, 10 /an. 
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Table 1.4. 

Effect of sample 
1120°C; flow rate - 

mass on reaction 
75 cc/min). 

kinetics (temperature - 

Run Number 

Sample Size 

(g) 

Carbon Content 
(wt%) 

Initial Rate 
(g CO/g ilmenite-s) 

1 

0.986 

20 

3.3 x 10- 4 

2 

0.498 

20 

3.1 x 10- 4 



Figure 1.13. Effect of carbon content on ilmenite reduction kinetics. 
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Another important finding is the influence of carbon type on the kinetics of the 
process. As shown in Figure 1.14, the reaction rate and even product distribution for 
two types of carbon are quite different. It seems that carbon forms with higher surface 
activity give a higher overall reduction rate. These preliminary results indicate that the 
gasification of carbon is an important (and possibly rate-determining) step in the 
overall process. Therefore, the use of a catalyst to promote this reaction should be 
investigated. One of the objectives of this study is to find a suitable catalyst to 
promote the gasification reaction. 

The chemical reactions taking place during the carbothermal reduction are as 
follows: 


C + FeTi0 3 - CO + Fe + Ti0 2 (1) 

CO + FeTi0 3 - C0 2 + Fe + Ti0 2 (2) 

C0 2 + C - 2CO (3) 

Reaction 1 is expected to start at the points of contact between carbon and ilmenite. 
However, the desirable path which will be emphasized in this project is a combination 
of Reactions 2 and 3. This path does not require solid contact and can use a wide 
variety of carbonaceous compounds for ilmenite reduction. 

Finally, a new temperature and concentration-programmed reaction is being 
developed to investigate the mechanism of the reduction reaction. This technique is 
suitable for studying the effect of various reaction conditions on rate using relatively 
few experiments. One example showing the effect of temperature on rate is shown in 
Figure 1.15. 
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Figure 1.14. Effect of carbon type on ilmenite reduction kinetics. 



Time (min) 

Figure 1.15. Temperature-programmed reaction results. 
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Abstract 

The use of chlorine to extract, reclaim, and purify metals has attractive possibilities 
for extraterrestrial processing of local planetary resources. While a complete cyclic 
process has been proposed for the recovery of metal lurgically significant metals and 
oxygen, this progress report focuses on the chlorination step of the cycle. An 
experimental apparatus for reacting refractory materials, such as ilmenite, in a 
microwave-induced plasma is being built. Complex equilibria calculations reveal that 
stable refractory materials can, under the influence of a plasma, undergo chlorination 
and yield oxygen as a by-product. These issues and the potential advantages for 
plasma processing in space are reviewed. Also presented in this report is a discussion 
of the complex equilibria program used in the analysis. 


Introduction 

On Earth, the use of chlorine to extract, reclaim, and purify metals has long been 
an attractive possibility. 1-5 In space, plasma chlorination has even greater possibilities. 
Those possibilities include the production of oxygen and chloride metal vapor. The 
latter can be either reduced or decomposed to its metal and chlorine states for 
recycling. In addition, a plasma chlorinator is versatile. It can be used to process a 
number of different ores, and it can also be used to decompose metal chloride vapor 
or waste products from other processes. 

Our immediate goal is to complete construction of the plasma reactor and begin 
testing. The unit will, for the next year, be the center of focus for our experimental 
program. The reactor has been designed to meet both our immediate needs and to 
serve as a service facility for other investigators participating in Center activities who 
may wish to determine if a plasma can enhance their research. 

In the same spirit of cooperation, we have assembled a software package for 
complex equilibria calculations and an extensive data base. The software package is 
based on Outokumpu's HSC program. The data base is oriented toward the extractive 
metallurgical industry and has already found considerable use in this work. The HSC 
program has been interfaced with LOTUS 123 for manipulation of data, Grapher and 
Energraphics for data plotting, and Inset to cut and insert results from any program into 
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another. The entire package is up and running, and we plan to put together a brief 
description regarding its operation. 

m 

Background 

The general advantages to chlorination processing are the large number of 
elements that form chlorides, the reactivity, volatility and solubility of metal chloride 
products, and the reducibility of metallurgically important chlorides to metal. 1 Almost 
all of the elements in the periodic table form chlorides, except for a few noble gases. 
In addition, chlorides tend to volatilize or have comparatively low melting points. Both 
physical phenomena provide a means for separation and subsequent refining. 

Metal chlorides, while easily produced by reaction of chlorine with the metal, are 
difficult to form from metal oxides. 2 This difficulty is usually overcome by employing 
a CO-CI 2 gas mixture or by reacting the oxide with Cl 2 in the presence of solid 
carbon. In carbochlorination, carbon monoxide acts as a reducing agent, forming C0 2 , 
while the chlorine oxidizes the metal, forming a volatile chloride. In this manner, a 
new surface is continually exposed for reaction. This process is certainly viable, and 
one worth considering for space application. However, there are technical complexities 
with carbochlorination that may be ignored on Earth, but cannot be overlooked with 
regard to space application. 

A number of chlorinated oxy-carbon compounds are stable at elevated 
temperatures. Landsberg et al. 6 conducted an extensive study of the relationship 
between various metal oxide carbochlorination reactions and C-CI-O compound by- 
products. They found and quantified 136 identifiable C-CI-0 by-product compounds 
produced during carbochlorination of oxides of titanium, zirconium, and aluminum at 
reaction temperatures between 600 and 1000°C. In an economic setting, such as that 
which can be anticipated for a space base, the loss of scarce reactants and oxygen in 
an unwanted by-product cannot be tolerated. 

"Cold” plasma chlorination can eliminate this problem, enhance reaction rates, 
minimize reactor size, and yield oxygen as a by-product. In addition, a plasma 
chlorination reactor is versatile. A single reactor could be used to process more than 
one ore, and it could also be used in the thermal reduction of chloride vapor. 

Approach 

"Cold" plasma processing of metallurgical materials is a natural progression of 
technological development within the extractive industry. High-temperature or 
"thermal" plasmas are used in many industrial applications in the extractive industry, 
including such areas as process heating and melting, thermal decomposition, smelting, 
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and refining. 78 Those applications rely primarily upon plasma heating to achieve 
conventional results that could be achieved with other heating sources. A "cold" 
plasma, unlike conventional heating sources, can significantly enhance reaction rates 
and form compounds, ions, and reactive species that cannot be produced by normal 
means. 

In spite of their name, "cold" plasmas are not necessarily cold. Depending upon 
the gas composition and pressure, a "cold" plasma can exist at molecular temperatures 
up to 2000 K. The term "cold" is attached to those plasmas wherein there is a 
significant difference between the molecular and electron temperatures, the latter often 
being thousands of degrees greater than the former. The reactive nature of a "cold" 
plasma is a direct result of the difference in the two temperatures. 

Both the reactive nature and the thermodynamic character of a "cold" plasma are 
represented by the sketch in Figure 1.16. In this sketch, the chemical potential of both 
reactants and products are plotted against an arbitrary variable to represent potential 
reactions. In a normal gas solid reaction, reactants combine to yield products and an 
overall reduction in the chemical potential of the system. In a "cold" plasma, energetic 
free electrons collide with gas molecules, leading to ionization and/or dissociation of 
those molecules. Depending on the initial gas composition, the excitation process can 
produce novel reactants, such as monatomic chlorine. 9 - 10 In addition, there occurs a 
substantial increase in the overall chemical potential of the reactants, which, in turn, 
allows a number of extra potential reactions, including the formation of metastable 
phases not attainable under non-plasma conditions. 

Reaction rates can be substantially increased in a "cold" plasma. The high 
chemical potential of ionized and dissociated molecules reflects the reactive nature of 
these molecules. The presence of broken bonds, partially filled orbitals, and 

unbalanced charge makes these molecules highly reactive and can substantially 
enhance kinetic processes. It is these characteristics of a "cold" plasma, and how they 
can be used to enhance extractive processes, which are the focus of this investigation. 

Results to Date 

Experimental Apparatus. The primary effort in the past eight months has involved 
design and construction of the experimental apparatus. A schematic diagram of the 
apparatus is presented in Figure 1.17. The system includes a microwave source, an 
applicator where the plasma is generated, a gas delivery system, a mass spectrometer 
for monitoring the extent of the reaction, an optical pyrometer for monitoring the 
temperature of the solid specimen in the plasma, a specimen holder that both rotates 
and allows for vertical translation of the solid specimen, and a vacuum system. 
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PLASMA 



Figure 1.16. Schematic representation of chemical potential for conventional and 
plasma-assisted reactions. 
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Figure 1.17. Experimental apparatus. Note: Mass spectrometer is rotated 90° and mass spec capillary 

sampling tube is not to scale. 
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The microwave generator was moved from a southside industrial site to the Arizona 
Material's Laboratory, where it has been set up and is now operational. The gas 
delivery ^system is in place; the specimen holder has been constructed, tested, and 
deemed operational. The reaction chamber, which is housed in the applicator, has 
been designed and built, while the applicator has been redesigned to accommodate the 
mass spectrometer. 

In addition to these activities, we are in the process of completing arrangements 
for specimen preparation. This activity has involved construction of a vacuum system 
for sealing specimens in silica ampules and setting up a sintering furnace. Specimens 
prepared from powders will be isostatically pressed and then sealed in a silica ampule 
under vacuum before sintering at approximately 1000°C. The solid specimens formed 
in this manner will be characterized using a BET apparatus to determine surface area 
and pore volume, with visual confirmation of those results provided through use of a 
scanning electron microscope. 

We are currently awaiting arrival of the mass spectrometer and mechanical pumps. 
The mass spectrometer, an EXTRELL C-50 Quadropole, was ordered last Fall. The 
instrument is currently being tested by the manufacturer, and the expected shipping 
date is March 17, 1989. Three corrosion-resistant pumps are on order and are 
expected to be shipped before the end of March 1989. Upon receipt of this equipment, 
we will complete construction of the apparatus and begin experimentation. 

Thermodynamic Calculations. In many systems, chlorination is limited by the 
dissociation of diatomic chlorine. In Figure 1.18, the equilibrium partial pressure of 
monatomic chlorine in equilibrium with the diatomic species is plotted assuming the 
sum of the partial pressures of both the monatomic and diatomic species equals 1 atm. 
Also in Figure 1.18, the partial pressure is compared to the partial pressure of 
monatomic chlorine that can be achieved in a "cold" plasma. The latter value, based 
on a total plasma pressure of 0.1 atm and total dissociation of diatomic chlorine, 
intersects the equilibrium-calculated curve at 1560 K. Depending upon the 
concentration of monatomic chlorine in the plasma, it is possible to enhance the 
chlorination kinetics greatly. 

The chemistry which takes place in a plasma is usually quite complex and involves 
a large number of elementary reactions which are initially generated by collisions 
between energetic free electrons and gas molecules. 10 - 11 While we cannot model the 
thermochemistry of the plasma at this early stage, we are in a position to estimate its 
potential impact on chlorination processing. The impact is estimated through complex 
equilibria calculations, using the computer program described earlier. 
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Figure 1.18. Comparison of the partial pressures of monatomic chlorine and hydrogen 
produced by thermal heating (total pressure 1 atm) and plasma (total 
pressure 0.1 atm). 

The chlorination of ilmenite has been modeled using the system presented in Table 
1.5. It has been assumed that, in plasma processing, diatomic chlorine is completely 
decomposed. Accordingly, the results presented below regarding the plasma represent 
the best-case scenario and do not take into consideration the effect of the plasma on 
other reactants and products. Moreover, it is important to remember that the results 
presented in this section are based on thermodynamic calculations and thus represent 
only what is energetically possible. Kinetic processes may also limit what can be 
achieved. 
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Table 1.5. System used for thermodynamic calculations. 


Vapor Phase 

Cl 2 and/or Cl, 0 2 , 

CIO, CI0 2 , FeCI, FeCI 2 , FeCI 3 , 

Fe 2 CI 6 , TiCI, TiCI 2 , Ti 2 Cl e , TiCI 4 , 

TiOCI, TiOCI 2 , 

TiO 


Solids 

FeTi0 3 

Ti0 2 

TiCI 2 

TiClg 

TiCI 4 

FeO 

Fe 2 O 3 

Fe 3 0 4 

FeCI 2 

FeCI 3 

FeOCt 



Ilmenite reacts with chlorine at 1100°C. The extent of this reaction is represented 
in Figure 1.19 as a function of the moles of chlorine available per mole of ilmenite. 
The results in the figure and the raw data reveal that the iron is volatilized as chlorides 
while elemental oxygen is released. The process continues until iron is eliminated, 
leaving only rutile (Ti0 2 ) remaining in the solid. At this point, the reaction stops for all 
practical purposes. 

The extent of the chlorination of ilmenite in a plasma can be enhanced, as the 
results shown in Figure 1.20 reveal. Even after all the iron has been volatilized, the 
chlorination process continues. Titanium chlorides are volatilized and elemental 
oxygen is released. A comparison of Figures 1.19 and 1.20 reveals that the chlorine 
may be better utilized in a plasma reactor. 

The potential for plasma chlorination is dramatically illustrated in Figure 1.21. 
Rutile, as noted earlier, is not readily chlorinated under standard conditions. The 
results in Figure 1.21 reveal that rutile, a highly refractory material, may readily be 
reacted in a plasma. It is this reactive nature that holds exceptional promise for the 
use of "cold" plasmas in the processing of local planetary resources. 

Summary and Future Work 

The investigation is in its initial stage. Considerable effort has been expanded in 
the design and construction of the plasma reactor. This initial work is expected to 
continue for three to four months. Upon completion of the apparatus, a number of 
preliminary experiments will be undertaken. 

The preliminary experiments will focus on the power absorbed by the plasma, 
which will be evaluated as a function of the gas content (chlorine versus inert gas 
content, as well as the choice of inert gas), the total pressure, reactor size, and power 
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Figure 1.19. Extent of reaction of ilmenite with chlorine under conventional heating at 
1100°C and total pressure of 0.1 atm. 



Figure 1.20. Monatomic chlorine in a plasma can readily react with ilmenite at 1100°C 
and a total pressure of 0.1 atm. 
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Figure 1 .21 . A plasma can effectively chlorinate refractory oxides. 


provided. The apparatus is equipped with meters to measure both the forward and 
reflected power, the difference being the absorbed power. The level of absorbed 
power reflects the extent of radical formation. The results of this work will be used to 
identify optimum conditions for the chlorination reaction. This work is expected to take 
no longer than two months. 

Chlorination studies will begin upon completion of the initial power studies. The 
kinetics of chlorination will be examined as a function of temperature, chlorine content 
in the gas, total pressure, power level, and plasma density. These experiments are 
expected to be in progress throughout the remainder of the year. 

Student Participants 

Mr. Daniel Bullard and Mr. Robert Ortega, both U.S. citizens, are working on this 
project. Daniel, who joined our research group in January, completed his 
undergraduate program at the University of Wisconsin at Milwaukee in December, 
where he received a double degree in Materials Engineering and Chemistry. Robert is 
in his last year of undergraduate study at the University of Arizona. Upon completion 
of his Bachelor's degree in May, he plans to commence graduate work. Robert began 
working on the project in December. 
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Previous work has shown that Fe° and 0 2 can be derived by electrolysis from 
silicate melt of a composition typical of lunar soils (Lindstrom and Haskin 1979). In the 
present study, our goal is to refine further the conditions necessary to optimize 
production and to determine efficiencies of production (how much product is derived 
for a given current) and purity of products. These depend on several factors, including 
potential imposed between electrodes, configuration and surface area of the electrodes, 
composition of the electrolyzed silicate melt, and oxygen fugacity. We have designed 
experiments to measure the dependence on these variables of three parameters that 
must be known before production by electrolysis can be optimized. These parameters 
are: 

Limiting Current - ffanode surface area (ASA), cathode surface area (CSA), 
distance between electrodes (d), electrolyte composition (X), potential 
between electrodes (E)] 

Actual Current - f[ASA, CSA, X, E, reduction potentials of electroactive species 
(E°), diffusion rates of electroactive species (D)] 

Efficiencies of Production - f(X, E, f0 2 ) • 

The limiting current is the maximum current the electrolyte can sustain at a given 
potential. It is related to the intrinsic conductivity of the bulk electrolyte. The 
conductivity has been measured previously for several compositions (Lewis 1985), and 
the limiting current can be calculated from these values. 

For a direct current to flow through the electrolyte, oxidation and reduction of 
electroactive species must take place at the anode and cathode. Therefore, the actual 
current depends on reduction potentials and concentrations of electroactive species in 
the melt, as well as the electrode surface areas, and in general will be somewhat less 
than the limiting current. Once this dependence is understood, the actual current can 
be optimized to the limiting current by adjusting ASA, CSA, and E. Because different 
reactions are occurring at the anode and cathode, the effect of ASA on the current is 
not necessarily the same as the effect of CSA. Therefore, it is necessary to measure 
their effects independently. We did this by using a three-electrode experimental set- 
up, one electrode to establish a reference potential and the other two as working 
electrodes. This permitted the total potential to be divided into the potential driving 
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the anodic reaction plus the potential driving the cathodic reaction. In this way, we 
have determined the dependence of current on ASA, CSA, E, and X. 

We have developed a method for measuring 0 2 yield and, therefore, 0 2 production 
efficiency. The amount of 0 2 produced is based on the difference in sample weight 
before and after the experiment. We have found a simple relationship which explains 
most of the variation in production efficiency for 0 2 . Mainly, efficiency decreases in 
proportion to the competing oxidation of Fe 2 * to Fe 3 *. 

Determination of efficiency of Fe production has proven more difficult because of 
the solubility of Fe° in the Pt wire loops used to hold the samples and the reduction of 
Si 4+ . However, we have estimated Fe efficiencies indirectly from the differences in 
current at various potentials. Experiments were designed to determine the difference 
between "background" current flow (related to reductions other than Fe) and total 
current flow. These values were determined by measuring current flow as a function of 
potential for both a silicate melt containing Fe and a silicate melt in which a molar- 
equivalent amount of electro- inactive MgO had been substituted for FeO. We found 
that Fe efficiency increases regularly with increasing FeO concentration and decreases 
with increasing potential. The decrease in efficiency with increasing potential is 
related to the increasing amount of Si 4+ that is reduced at higher potentials, suggesting 
that the purity of the Fe product also decreases with increasing potential. 

In summary, we have determined the three relationships above for several 
compositions over a range of potentials from 0 to about 1 volt. With this information, 
the three equations can be solved simultaneously to optimize ASA, CSA, E, X, and f0 2 
for optimum production rate or production efficiency. Furthermore, because these 
relationships are general, rather than specific to a preconceived electrolytic set-up, 
these data will be applicable to whatever electrolytic set-up further research finds to 
be most practical. 

An abstract describing these results which was published in Lunar and Planetary 
Science XX follows. 

Some Effects of Composition on Efficiencies for 
Production of 0 2 and FeP From Silicate Melts by 
Electrolysis , R. 0. Colson and L. A. Haskin, 

Department of Earth and Planetary Sciences and 
McDonnell Center for the Space Sciences, Washington 
University, St. Louis, Missouri 

Several experiments have been conducted aimed at determining the efficiencies for 
production of 0 2 and Fe metal by electrolysis of silicate melts. Efficiency of 0 2 
production [defined as the moles of 0 2 produced/(moles of electrons passed through 
' the melt/4)] is calculated from the amount of 0 2 produced and the measured total 
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current flow. The amont of 0 2 produced is determined from the difference in sample 
weight before and after the experiment. Efficiency of Fe metal production was 
determined by a combination of theoretical treatment and experiment. Experiments 
were designed to determine the difference between "background” current flow (related 
to reductions other than Fe) and total current flow. These values were determined by 
measuring current flow as a function of potential for both a silicate melt containing Fe 
and a silicate melt in which a molar-equivalent amount of electro-inactive MgO has 
been substituted for FeO. Because the duration of these Fe-efficiency experiments is 
short and the species being reduced are not significantly depleted, these estimates 
correspond to the case in which the melt being electrolyzed is continuously "refreshed" 
and well stirred (except that oxidation products from the anode are not reduced at the 
cathode). 

Efficiencies of production will be influenced by the presence of reducible or 
oxidizable species other than the species of interest. For example, 0 2 production is 
expected to be made less efficient by the concurrent oxidations of Fe 2 * * Fe 3 * and 
Fe° .* Fe 2 * (we have found a significant amount of Fe° to be soluble in the silicate 
melt, similar to solubilities of Ni°, Co°, and Zn° found by Semkow et al. 1 ). Fe° 
production is expected to be made less efficient by reductions such as Fe 3 * ■* Fe 2 *, 
Si( 4 +) -* Si°, Cr 3 * *♦ Cr 2 *, and Ti 4 + •* Ti 3 +; however, Cr and Ti are not present in the 
experiments of this study. We expect the efficiencies of 0 2 or Fe production to 
decrease as the relative amounts of competing species increase and as electrode 
potentials approach the reduction potentials of the competing reactions. Therefore, it is 
useful to look at efficiences as a function of (1) electrode potential, (2) oxygen fugacity, 
and (3) melt composition (in specific, Fe concentration). 

0 2 Efficiency.— 0 2 efficiencies as a function of electrode potential are shown in 
Figure 1.22. The low point in the oxygen efficiency curve is near the reduction 
potential for Fe 3 * which, as expected, is the region of minimum efficiency. (Based on 
Fe 3 * reduction potentials, we estimate about 35-65% of the Fe 2 * is oxidized at these 
potentials.) At low potentials, a small amount of Fe 2 * is oxidized to Fe 3 * (we estimate 
5-10%) and at higher potentials the amount of Fe 2 * that is oxidized is small compared 
to the amount of oxygen because oxygen is much more abundant in the melt than Fe. 

0 2 efficiency as a function of f0 2 at constant potential and constant total Fe 
concentration is shown in Figure 1.23. The efficiency is lower at lower f0 2 , consistent 
with the higher concentration of Fe 2 * at the lower f0 2 and, therefore, a larger fraction 
of current is related to the oxidation of Fe 2 *. For this change in oxygen fugacity, the 
fraction of current attributable to Fe 2 * - Fe 3 * increases by 65%, and calculations based 
on the report of Sack et al. 2 indicate an increase of 55% in Fe 2 * concentration. If we 
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assume that all the current is related to oxidation of either oxygen ions or Fe 2 *, we can 
choose the equation %Fe 2+ /%0 2 - X(Fe 2+ )/X(02-) to represent relative amounts of the 
current attributable to each oxidation, where %Fe2* is the percentage of the current 
attributed to the oxidation of Fe 2 * and X(C) 2 ') is the effective concentration of 
oxidizable oxygen (which we assume is constant for the series of compositions between 
diopsidic and hedenburgitic melts used in this study). The equation can be modified 
(assuming %Fe 2+ + %0 2 • 100%) to yield %O 2 /(100 - %0 2 ) « X(0 2 ~ )/X(Fe 2 * ). This 
relationship predicts that as Fe 2+ concentrations increase, 0 2 efficiency will decrease. 
This is the relationship observed and shown in Figure 1.24. 

Fe Efficiency.— Figure 1.25 plots cathodic current as a function of electrode 
potential for several compositions. The "background" current is for a composition in 
which MgO has been substituted for FeO. This current is related primarily to the 
reduction of Si 4 *. If we assume that no current is lost to electronic conduction, 3 then 
the figure gives a rough estimate of the relative amounts of the total current that can be 
attributed to the Fe 2+ and Si 4+ reductions at various potentials and Fe concentrations. 
However, for these compositions and potentials, most, if not all, of the reduced Fe is 
dissolved in the melt or Pt electrodes and does not exist as (readily extracted) pure 
metallic Fe. Threshold potentials for the precipitation of pure metallic Fe can be 
calculated from the data presented in Grove. 4 Threshold potentials so calculated are 
plotted as a function of Fe concentration in Figure 1.26. Comparison of Figures 1.25 
and 1.26 indicates that in the range of potentials we have studied, pure metallic Fe can 
only be precipitated at concentrations of Fe > 35 mole%FeO and, at those 
concentrations, the amount of Fe reduced will be very large compared to the amount of 
Si 4 * reduced. These results indicate that high efficiencies of 0 2 and Fe production can 
be achieved by electrolysis, but not in bulk melt of the same composition. High 
efficiencies of 0 2 production are favored by low Fe concentrations, and high Fe° 
production efficiencies are favored by high Fe concentrations. 

Acknowledgments.— We gratefully acknowledge partial support for this work from 
UA/NASA Center for Utilization of Local Planetary Resources. 
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During the time period from September 1, 1988 to February 28, 1989, work has 
progressed in the following areas. Lunar simulant material has been fabricated to 
simulate Apollo 11, Apollo 15, and Apollo 16 soils. Table 2.1 shows the approximate 
compositions which these soils simulate. Samples of each soil simulate were fabricated 
into right circular cylinders of approximately 1.27 cm diameter by 1.27 cm length. 
Each sample was then placed into a microwave reaction cavity and heated using 2.45 
GHz radiation to a particular processing temperature at which it was held for a 
prescribed length of time. Figures 2.1 and 2.2 show examples of two microwave 
heating programs for these simulant samples. Upon completion of thermal processing, 
each sample had its top and bottom surfaces polished to a surface roughness of 
approximately 4 \l inches RMS. 


Table 2.1. Simulated compositions (weight percent). 



All 

(10089, 4) 

A15 

(15999, 126) 

A16 

(69999, 75) 

Si0 2 

42.03 

46.52 

45.26 

Ti0 2 

7.48 

1.40 

0.58 

ai 2 o 3 

13.59 

17.10 

26.22 

FeO 

15.74 

12.05 

5.82 

MnO 

0.20 

0.17 

0.07 

MgO 

7.86 

10.42 

6.39 

CaO 

11.98 

11.17 

14.76 

Na 2 0 

0.44 

0.43 

0.45 

K 2 0 

0.14 

0.20 

0.13 

Cr 2 0 3 

0.30 

0.27 

0.12 

p 2 o 5 

0.11 

0.20 

0.13 

s 

0.13 

0.07 

0.07 


Samples were then thermally shocked using an apparatus consisting of a tube 
furnace vertically oriented so that a sample could be suspended in it while being 
heated and then be allowed to drop, after some prescribed- amount of time, through the 
furnace into a water quench. 
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Figure 2.1. Time versus temperature for samples 2D and 4D. 
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Figure 2.2. Time versus temperture for samples 9E-1 and 9E-2. 
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Initial work has been reported by Meek et al. (1988) on thermal shock properties of 
microwave sintered lunar simulant materials. The mechanical testing of these samples 
was carried out using conventional ASTM (American Society for Testing Materials) 
techniques. There is some question, however, about the procedure, so it was decided 
to look again at samples using a different technique, one which employs hardness data 
to determine the fracture toughness of a material and then relates this to thermal shock 
resistance. Some advantages of this technique are its low cost for analysis, a sample 
can be tested in many locations, and thermal shock resistance as a function of initial 
sample orientation in the microwave field can be determined. 

Currently, we have prepared samples for testing at the Y-12 facility of Oak Ridge 
(operated by Martin Marietta) and are awaiting results. At the same time, other samples 
are currently being investigated using scanning electron microscopy to determine their 
microstructures. 

An abstract for a paper on this work has been submitted for presentation at the 9th 
Princeton Conference on Space Manufacturing. The only uncertainty concerning this 
paper is whether we will receive data from the Y-12 facility in time to present it in this 
paper. If we do not, the paper will focus mainly on microstructural analysis of the 
samples as a function of processing technique. 

For the next 16-month period, we propose to continue the work on lunar simulants, 
begin work on actual lunar sample material (Apollo 11, 15, 16), and begin work on 
carbonaceous chondrite material. Samples will be prepared in the same manner as the 
simulants. They will be thermally processed by conventional means and by using 2.45 
GHz radiation. Their mechanical properties and thermal shock properties will be 
determined and then related to their observed microstructure. 

The proposed budget for March 1, 1989 through June 31, 1990 is as follows: 

I. Salaries 



A. 

Graduate Student 

$15,872 


B. 

Advisor 

5,000 


C. 

Shop 

1,000 


D. 

Secretarial 

2,372 


Benefits at 28%% of I.B-D 

1,820 


Subtotal 

$24,192 

II. 

Materials 

1,000 

III. 

Travel 

800 


Subtotal 

IV. Overhead at 38.5%% of I, II, and III 


25,992 

10,008 


TOTAL 


$36,000 
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Abundances of Volatile-Bearing Phases in Carbonaceous Chondrites and Cooling ^ £ 
Rates of Meteorites Based on Cation Ordering of Orthopyroxenes 

J. Ganguly 

Department of Geosciences 
The University of Arizona -- 

The "deliverables" for the grant received during this period are (1) a 
comprehensive proposal for work during the period 3/1/89 to 2/28/90 and (2) results of 
preliminary calculations of volatile abundances in carbonaceous chondrites. 

Both objectives have been met. The results of the calculations of the volatile 
abundances will be presented at the Ninth SSI Princeton Conference on Space 
Manufacturing, May 11, 1989. In addition, the method (Ganguly 1982) has been refined 
for the calculation of cooling rate on the basis of cation ordering in orthopyroxenes, 
and it has been applied to the derivation of cooling rates of some stony meteorites. 
Evaluation of cooling rate is important to the analysis of condensation, accretion, and 
post-accretionary metamorphic histories of meteorites. The method of "orthopyroxene 
speedometry" is widely applicable to meteorites and would be very useful in the 
understanding of the evolutionary histories of carbonaceous chondrites, especially since 
the conventional metal lographic and fission track methods yield widely different results 
in many cases. The results will be presented at the Lunar and Planetary Science 
Conference XX, 1989. The following abstracts summarize the major conclusions of the 
volatile abundance and cooling rate calculations. 

Theoretical Predictions of the Nature and Abundances 
of Volatile Bearing Phases in Carbonaceous Chondrites, 

Jibamitra Ganguly, Department of Geosciences, 

University of Arizona, and Surendra K. Saxena, 

Department of Geology, Brooklyn College 

The overwhelming majority of asteroidal bodies consist of carbonaceous 
chondrites. Thus, the study of carbonaceous chondrites is very important to the 
assessment of resource potentials of near-earth asteroids. The carbonaceous chondrites 
are important potential resources for volatiles which may be used for fuels and 
propellants. The volatiles are stored in the minerals as structurally bound components 
in the organics, and in the possible cometary cores. The hydrous minerals have been 
characterized as a mixture of phyllosilicates, but further characterizations by 
conventional methods have proved to be extremely difficult due to the extremely fine- 
grained nature of the material. However, these phases need to be properly 
characterized if we are to efficiently release their volatiles for practical uses. 
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In this work, we have calculated at several P-T conditions the nature and 
abundance of the volatile bearing phases by minimizing the Gibbs free energy (G) of 
the system Mg-Fe-AI-Si-C-H 2 -C> 2 > constrained by the bulk composition of the Orgueil 
Cl chondrite. The minimization was carried out using the method of Lagrangian 
multiplier and a thermochemical data base consisting of about 100 minerals and 70 
vapor speices. Further, the effects of solid solutions (ideal and nonideal) have been 
taken into account. Due to the lack of adequate data, we have assumed that the 
phyllosilicates behave as ideal solutions. 

The results show that a minimum pressure of about 1 bar at 500 K is needed to 
form the hydrous phases and that the primary hydrous minerals are antigorite and talc. 
Some amount of chlinochlore (~ 2 mol%) may also be present. The amount of water 
structurally bound in the minerals should be about 5-10 wt% of the condensed 
inorganics. The amount of crystalline graphite may be very small, within a few mol%. 
There is also abundant magnetite. The relatively high pressure required to form the 
hydrous phases could preclude their formation by direct condensation from solar gas. 

The vapor phase in equilibrium with these crystalline phases consists primarily of 
H 2 0 (~ 70 mol%), CH 4 (~ 20 mol%), and C0 2 (~ 10 mol%). This may be an 
approximate composition of the possible cometary core if we assume that the core 
represents the trapped volatiles, the composition of which was controlled by reaction 
with the minerals and other condensed phases. 

There are considerable uncertainties in the thermochemical data of the 
phyllosilicates, especially of the Fe-end members and in the thermodynamic mixing 
properties of their Fe-Mg joins. Thus, the above results must be treated with caution. 
Further studies on the refinement and systematization of the thermochemical properties 
are in progress. Also, we need to take into account sulfur and organics in free energy 
minimization calculations to develop a more reliable model of the mineralogical 
characteristics of carbonaceous chondrites. This work is to be published in the 
Proceedings of the May 1989 SSI Conference (Appendix E). 

Acknowledgment'. NASA Grant NAGW-1332. 

Fe^ + -Mg Ordering in Orthopyroxenes and the Cooling 
Rates of Meteorites , Jibamitra Ganguly and Kunal 
Bose, Department of Geosciences, University of 
Arizona, and Subrata Ghose, Department of Geological 
Sciences, University of Washington 

The cooling rates of meteorites have been a subject of central importance in the 
understanding of their origin and evolution. However, the commonly used methods for 
the determination of cooling rates of meteorites often yield conflicting results. A 
potentially powerful recorder of the cooling rates is the state of Fe-Mg ordering in 
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orthopyroxenes (OPx), which can be found in any stony and stony-iron meteorite, and 
as silicate inclusion in some iron meteorites. However, despite its many applications in 
terrestrial problems, this method has not yet been used to derive quantitative 
constraints on the cooling rates of meteorites. 

The "cation-ordering speedometry" involves reconstruction of the ordering state of 
a sample at its peak temperature and numerical simulations of the change of site 
occupancies as a function of cooling rate (Ganguly 1982). A substantial amount of 
thermodynamic and kinetic data, which are required for these calculations, is available 
in the literature (Ganguly 1982, Anovitz, et al. 1988). An example of the simulated 
evolution of M2 site occupancy is shown in Figure 2.3, using an OPx from the Steinbach 
meteorite for which the measured XFe(M2) - 0.304 [XFe - Fe 2+ /(Fe 2+ + Mg)]. The 
cooling is assumed to follow the form 1/T - 1/T 0 + nt, where T 0 is the peak 
temperature and n is a cooling time constant. The observed site occupancy is exactly 
reproduced as the quenched ordering state in a simulated evolution, which is 
characterized by T-t path with n - 1.5(10- 12 ) yr 1 . Similar calculations are carried out 
for orthopyroxenes from Shaw and Johnstown meteorites (Table 2.2). The cooling 
rates, given by nT 2 , are strictly applicable around the closure temperatures (T c ) and 
are quite insensitive to errors in the estimation of T 0 . Any other cooling law can also 
be used. 

Figure 2.4 shows the effect of cooling rates on quenched M2 site occupancy for 
three bulk XFe values. For comparison, all XFe(M2) values are relative to those for a 
cooling rate of 150° C/my for the different bulk compositions. It is clear that the 
retrieved cooling rates are very sensitive to errors in the determination of site 
occupancies of natural samples, especially for Fe-poor compositions. A precision of 
0.002 in the site occupancy determination is a difficult but attainable goal in single 
crystal X-ray analysis. This will limit the precision of the retrieved cooling rates of 
meteoritic orthopyroxenes to a factor of 5-10, depending on the composition and 
cooling rate. 

The site occupancies of the meteoritic orthopyroxenes used in this work lack the 
high precisions required for reliable cooling rate calculations. Nonetheless, it is 
important to note that, despite their uncertainties, the inferred cooling rates for Shaw 
and Steinbach orthopyroxenes are in good agreement with those derived from 
metal lographic properties of these meteorites (Taylor et al. 1979, Reid et al. 1974). In 
contrast, the 24 Pu track records were used to derive (Pel las and Storzer 1977, Pel I as et 
al. 1978) a cooling rate for the Shaw meteorite below 350°C, which is three orders of 
magnitude slower than that shown in Table 2.2. Thus, either there was a drastic change 
of cooling rate of the Shaw meteorite between 400 and 350° C, or the 24 Pu result 
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Figure 2.3 


Table 2.2. Summary of cooling rate calculations of meteoritic orthopyroxenes. 


Sample 

Fe 2 * /(Fe 2 *+Mg) 
Total 

M2 

O 

O 

n.yr 1 

(dT/dt) Tc , 

°C/my 

Site Occupancy 
Measurement 

Johnstown 

0.175 

0.324 

400 

4E-9 

1.8E3 

X-ray 

(Dodd et al. 1975) 

Steinbach 

0.161 

0.304 

304 

1.5E-12 

0.5 

X-ray 

(Reid et al. 74) 

Johnstown 

0.237 

0.451 

318 

1.5E-12 

0.5 

X-ray a 


a Preliminary results from this work. 
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requires a careful reevaluation. Although the Johnstown meteorite is shocked (Mori and 
Takeda 1981), the high degree of Fe-Mg order in OPx suggests that the magnitude of 
shock was <450 Kb (Dundon and Hafner 1971). These results illustrate the potential 
applications and importance of the "cation ordering speedometry," which should be 
most widely applicable in the study of cooling rates and shock histories of meteorites 
owing to the widespread occurrence of orthopyroxenes. 

Acknowledgments'. NSF Grant EAR 8206656 and NASA Grant N AG W- 1332. 
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The first six months were devoted to (1) reviewing the chemical and physical 
characteristics of ilmenite in lunar material and (2) developing high-precision, low- 
abundance analyses for platinoids and other possibly important elements in space 
materials using an induced coupled plasma mass spectrometer. 

Lunar Ilmenite 

There is considerable literature on the chemistry of lunar ilmenite and the 
differences between lunar and terrestrial ilmenite. It is clear that terrestrial ilmenite is 
a poor analog to lunar ilmenite and should be avoided as a substitute material in any 
experiment since its physical and chemical properties will be different. If substitute 
material is needed for lunar ilmenite, it would be preferable to use synthetic ilmenite. 
A problem that was discovered during the literature review is that lunar regolith, which 
would be the starting material for ilmenite beneficiation, is not that well characterized. 
It is critical to know if ilmenite is free, the size of average ilmenite grain, and the size 
of average regolith material in order to be able to physically separate the ilmenite from 
the regolith. Consequently, we are in the process of asking NASA for lunar regolith so 
that we can characterize the material with ilmenite beneficiation in mind. 

High-Precision Analyses of Platinoids 

The techniques have been developed to analyze solutions which have a few parts 
per billion platinoids with reproducibility that is better than 7%. This will allow us to 
characterize any planetary material to determine if platinoids could be "mined" as a 
subproduct. 

Future Work 

In the next year we will characterize the physical characteristics of lunar regolith 
and continue developing the techniques to analyze small amounts of planetary material 
for economically viable elements. 
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This work involves the theoretical calculations of rocket performance of "exotic" 
fuels at various operating conditions such as chamber pressure, pressure ratios, and 
oxidizer-to-fuel ratios. By "exotic" fuels, we mean using materials that may not 
normally be used as fuels here on Earth due to low performance. The majority of the 
work is currently being done on a VAX using the CET 86 program by Gordon and 
McBride. The CET 86 program will be installed on our IBM PS/2 Model 80 to reduce 
computer costs. We have also received the JPL Solid Propellant Theoretical 

Performance program for the IBM, which handles solid fuels better than the CET 86 
program. This program was developed by the Air Force and has been made available 
to us for restricted use. In addition, we also have received the AVCO version of the 
program. These enable easy comparison. 

With the high cost of putting payloads into orbit, the ability to reduce the size and 
weight of the payload is highly beneficial. Since, in any space vehicle, the largest part 
of the vehicle consists of propellant, it is hoped that by finding space materials suitable 
for use as propellants, the overall size and weight of the vehicle can be reduced. To 
this end, this work is attempting to find propellants that will give acceptable 
performance, yet be producible or readily available in space scenarios. Another 
possibility is the use of materials, such kevlar and rubbers, that would be taken for 
other purposes. Also considered is the need for the propellant to be easily stored over 
long periods of time while ideally avoiding the use of heavy refrigeration systems. 

Although extensive research has been done on Earth-transported propellants that 
have high performance characteristics, very little work has been done on propellants 
that give less than ideal performances. In additon to studying such overlooked 

propellants, we need to consider low chamber pressures, to help reduce weight, and 
high pressure ratios which can be achieved in atmospheres such as on Mars or the 
Moon. Some propellants considered unacceptable for use on Earth could be ideal for 
use elsewhere since trade-offs of performance for availability, storability, and cost are 
acceptable within the scope of this research. 

We are approaching this problem by setting up the specific programs first on a 
VAX and on a PC. For the propellants I am presently considering, I am varying the 
chamber pressure and the oxidizer-to-fuel ratios for a variety of pressure ratios. Since 
just the present list of propellants will generate over 22,000 sets of data, we plan to 
present the data graphically, as well as in tabular form. After all the propellants have 
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been analyzed, we will be able to determine which ones (at which conditions) best 
meet our requirements. 

We have set up the CET 86 program by Gordon and McBride on the VAX and have 
obtained results for a number of the propellant combinations being analyzed. A sample 
of these results is included (Appendix B). Several programs have been written to 
reduce the large amounts of data generated to include only the information needed, 
such as specific impulse, temperature, etc. We are currently working on a program to 
graphically display some of these data. Most of the initial runs were "verification" runs 
on conventional propellants. Nevertheless, one "exotic" extraterrestrial propellant is 
also shown in this appendix. 

This work involves analyzing propellants at various operating conditions and 
determining which will give acceptable performance and yet be easily storable for long 
periods of time and be readily available in space missions. Work still to be done 
consists of completing the analysis of the remaining propellants and organizing the 
results into a more useful format. To assist in organizing the data, a commerical 
plotting program, such as the MathCAD program, would yield results in a more timely 
manner and make the development of a plotting program unnecessary. We will procure 
this and install it in 1989-90. 

Jennifer Kares, an undergraduate student, has been assisting in the analysis of the 
fuels. Jeff Kahl, also an undergraduate student, has been assisting in the development 
of the various programs needed to process the data. These two are volunteering their 
time for academic credit, but at no cost to NASA at this time. 
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Abstract 

The development of construction materials such as concrete from lunar soils 
without the use of water requires a different methodology than that used for 
conventional terrestrial concrete. In this research project, a unique approach will be 
attempted that utilizes factors such as initial vacuum and then cyclic loading to enhance 
the mechanical properties of dry materials similar to those available on the moon. The 
application of such factors is expected to allow reorientation, and coming together, of 
particles of the materials toward the maximum theoretical density. If such a density 
can provide deformation and strength properties for even a limited type of construction, 
the approach can have significant application potential, although other factors such as 
heat and chemicals may be needed for specific construction objectives. 


Introduction 

Significant progress has been made in the area of space exploration in the past 
three decades. However, little has been done for development of construction 
materials for long-term shelters and work stations for man on other planets. For large- 
scale productivity on other planets, buildings and other structures are indispensable. 
Because of payload limitations on spacecrafts and the high cost of transportation, a 
prudent and efficient approach to construction on other planets requires that every 
effort be made to utilize the materials available on the planets themselves. The 
technique that is being developed in the Department of Civil Engineering and 
Engineering Mechanics is intended to utilize extraterrestrial "soils" to develop 
construction materials for building safe and serviceable structures. 


Approach 

Samples of lunar soil simulants of specimen size (3.0 x 6.0 inches) will be mixed in 
different proportions and will be tested under various levels of initial density and 
vacuum. Compacting the soil at different magnitudes of vacuum will give an indication 
of the effects of vacuum on lunar soil compactability. Then, cyclic loading, causing 
compression of the initial specimen, will be applied in a specially designed chamber 
using an MTS test frame. The experimental results will be plotted in terms of relations 
between initial density, initial vacuum, and number of loading cycles. 
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The material characteristics of the compacted samples (stress-strain response, 
angle of internal friction, and "apparent” cohesion) will be evaluated to assess the 
deformation and strength capabilities of the resulting compacted specimens. 

Indeed, application of only compression cycles may not be able to provide the 
required deformation-strength characteristics for (all) construction objectives, and 
additional research to include other factors such as use of heat and chemicals would 
then be considered in subsequent research. This initial research would, however, 
provide the limit of compaction that can be achieved with cycles of compression 
without the use of water. 

Results to Date 

Samples of lunar soil simulant in grain sizes from 2 millimeters to less than 75 
microns have been produced by crushing terrestrial basalt rocks. Approximately 0.25 
cubic foot of lunar soil simulant, 30 percent of which is ultra-refined, is available for 
testing purposes. 

A cylindrical vacuum chamber 15.0 inches high by 3 inches inside diameter has 
been machined from a single block of aluminum; design and fabrication of this device 
are supported by a grant from the Small Grants Program of the University of Arizona. 
The chamber is equipped with a ram-piston assembly to provide a means of loading 
the sample and is lined with Teflon to prevent cold-welding of the soil grains to the 
aluminum. Two ports in the chamber allow application of vacuum for both the top and 
bottom of the soil sample. Pressure is measured by a digital vacuum gauge with two 
Pirani PG-3 sensors attached to the chamber. The system is attached to a mechanical 
roughing pump and a diffusion-type vaccum pump capable of reducing the chamber 
pressure to 10 torr. Cyclic load is applied using an MTS hydraulic testing device 
which allows variation in load, cycle duration, and number of cycles and can 
automatically measure the displacement of the sample during loading. 

Planned Research During the Next Year 

A pilot series of tests will first be performed to establish the maximum compaction 
that can be achieved under different initial vacuums and density for a specimen of 
given initial height. 

Then, a series of tests will be performed with three to four values of initial vacuum 
(V 0 ) and density (p 0 ). For each, the specimen will be subjected to cycles of 
compression loading. The results will be plotted in terms of variation of density with 
number of load cycles (N) for each Vg and pg. These results will be used to develop 
nondimensional relations between changing density as affected by V 0 , p 0 , and N. 
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A modified test device will be developed so that the compacted specimen can be 
tested with increasing axial stress, cr , , (near) unconfined conditions, i.e., zero confining 
pressure o r 3 and different values of ff 3 . The results will provide stress-strain relations* 
in terms of cr-, - cr 3 vs. e, (axial stress). The deformation moduli and the strength 
parameters (angle of friction <f> and cohesion C) can then be found from these stress- 
strain curves. 

The above information will allow evaluation of the deformation-strength 
characteristics of the compacted specimens and the range of loading due to 
construction that such material is able to withstand. 

It is felt that the above approach may not lead to materials with strengths required 
for (high) loads that may arise in common construction. Hence, the second phase of 
research would involve the use of factors such as heat and chemicals in addition to the 
cyclic compression loading. Acquisition of a furnace for this purpose is being 
considered at this time. 
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This program was undertaken in order to understand the complex nature of 
interdependent performance in space missions. As the first step in a planned sequence 
of progress, a "spread sheet" program has been developed to evaluate different 
fuel/oxidizer combinations for a specific Martian mission. This program is to be linked 
with output attained using sophisticated software produced by Gordon and McBride. 
The data must be input by hand at this time. The programming to date makes use of 
11 independent parameters. While the work is somewhat developed, some 

modifications such as the inclusion of aerobraking, in-situ propellant, output 
presentation, and automation of the data feed are still necessary. The program should 
prove to be quite valuable, even under its early stages of development. 

Optimization is essential when faced with the incredible magnitude of costs, risks, 
and benefits involved with space exploration. We must decide not only what to use, 
but how to utilize it most effectively. Making the most of what we have and what we 
find is the essence of our NASA Center. A system of weights needs to be devised on 
which to measure our options. It is the goal of this work to devise a Figure of Merit 
(FoM) on which different choices can be presented and made (Fig. 3.1). The plan is to 
model typical missions to Mars, identify the parameters, and vary them until the best 
one is found. Initially, most of the focus is placed on propellant selection. 

There has been some work done in the past and some currently under investigation 
which will prove crucial to our goals. Dr. Frisbee at JPL has been working on 
autonomous oxygen and carbon monoxide production from the Martian atmosphere for a 
number of years. He has provided his most current papers on the subject. They 
include the modeling information for power requirements, production rates, and system 
mass breakdown. There is a group willing to provide us with very reliable incremental 
velocity requirements for many different scenarios. There is a group in Austin, Texas, 
who have been doing feasibility studies for different lunar missions. They have 
neglected some important parameters regarding propellant selection, but their input still 
may prove useful. The tasks ahead are open for investigation. 

The study utilizes current spread sheet technology and effectively demonstrates 
how changing a single parameter influences the feasibility of a mission. An IBM PC is 
currently being used to run LOTUS 123 Version 201. The plan is to keep the 
programming simple at first, building on the complexity with time. The first run of the 
program is for a specific mission, starting at LEO, making a landing on Mars, and from 
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Figure 3.1. "Figure of merit" optimizations. 
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there placing a given payload into a Martian satellite orbit. The rocket staging was 
developed from the velocity requirements using Newton’s laws and a minimum-energy 
transfer. This approach will utilize CET 86, Gordon and McBride’s software for 
unconventional fuel performance. The current input parameters are made up of those 
placed on the fuel/oxidizer combination, the nozzle, the structural mass coefficient 
(tanks, pumps, refrigeration units, motors, etc.), and the payload itself. To date, they 
include the selection of the fuel and the oxidizer, whether the expansion is modeled as 
frozen or equilibrium flow, the mixture ratio, the chamber pressure, the ambient 
pressure, the expansion ratio, the nozzle’s material density, its average thickness, the 
structural coefficient, and the payload mass. From this information, some useful 
quantities are found for each stage of the rocket and for the vehicle as a whole. They 
are made up of the following: exhaust velocities, mass flow rates, vehicular thrusts, 

throat areas, exit plane areas, nozzle coefficients, payload coefficients, initial masses, 
final masses, structural masses, fuel/oxidizer masses, and FoM. Currently, the FoM for 
each stage, as well as for the vehicle, is defined as the total initial mass divided by the 
mission payload. This will be modified to include risk, reliability, support hardware, 
etc. The real beauty of this approach to programming is the ease of immediate 
updating of the results. 

To date, some runs have been accomplished. The last two months have been 
spent developing the basic logic for the program. A few additional parameters must be 
examined before the more visual aspects of the work can be effectively presented. 
The first is automating the data entry. All changes in the input have been entered by 
hand; a more efficient manner is under investigation. The LOTUS software is equipped 
with "macros" which allow complex, yet tedious tasks to be performed. The results 
should prove to be very powerful, especially when considering that over ten thousand 
points of propellant data have already been generated using CET 86. Having a 
subroutine to enter the data into the main program is not only feasible but also 
necessary for computer processing. The next step is to include an "aeroshell" in the 
model. The huge mass savings will prove to be a major factor in the minimization of 
the FoM. It will also make the mission profile much more realistic. Once these tasks 
have been accomplished, the real motivation for this initial spread sheet can be more 
effectively explored. In-situ propellant production will also start off simply, with the 
extraction of oxygen from the Martian atmosphere. These additions will make the data 
more presentable, both to those interested in finding an optimal solution and to others 
interested in viewing families of curves generated from other output data. 

The research and complexity of the program are developing at a steady rate. The 
first major run of the spread sheet will come after the inclusion of "aerocapture" into 
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the model. The basic program will continue to evolve. The FoM must also contain a 
reliability factor. Leakage rates for all Earth-carried oxygen need to be included. The 
structural factor needs to be calculated using fuel density, chamber pressure, and 
cryogenic requirements. The different possibilities for in-situ propellant production 
must be examined. Oxygen can be produced from the atmosphere or from the polar 
ice cap. Other possibilities include hydrocarbons from water and C0 2 . Further 
complexities will enhance the value of the work but will slow the computation speed. I 
have been in communication with the people at LOTUS regarding their Version 3; it is 
due to be released in the second quarter of 1989. Their newsletter claims that it meets 
their goals for speed and performance and, "123 Release 3 features significant power 
and performance enhancements including true three-dimensional spreadsheets, external 
database access, new recalculation techniques and graphics and printing 
improvements." We should be able to purchase an upgrade of the current version for 
about $200 and run the software on either of the SUN386i computers currently at the 
Center. The three-dimensional version would be a very useful means of storing both 
inputs and outputs. 
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Abstract 

Three specific projects were initiated with late starting dates. The first and 
second projects, "Automatic Control of a Propellant Mixer" and "Computer Simulation 
of Engineering Processes," were started January 15, 1989, and the third project, 
"Control for Chaotic Mixing," was started February 1, 1989. 

The first project is concerned with the automation of a small, solid-propellant 
mixer located on campus in the AME labs building. Temperature control is under 
investigation. A numerical simulation of the system is under development and will be 
tested using different control options. Control system hardware is currently being put 
into place. 

The second project deals with the construction of mathematical models and 
simulation techniques for understanding various engineering processes. Computer 

graphic packages will be utilized for better visualization of the simulation results. 

The third project is concerned with understanding and, hopefully, improving the 
mechanical mixing of propellants. Simulation of the mixing process is being done in 
conjunction with the second project. We will use this simulation capability to study 
how one can control for chaotic behavior to meet specified mixing requirements. An 
experimental mixing chamber is also being built. It will allow visual tracking of 
particles under mixing. The experimental unit will be used to test ideas from chaos 
theory, as well as to verify simulation results. This third project in the basic research 
category has clear applications to extraterrestrial propellant quality and reliability. 


Introduction 

In the production of solid propellants, components must first be mixed under 
specified conditions in order to produce a uniform product. To date, the batch 
processing of propellants is still largely a process with a human operator in the control 
loop. Not only is this unacceptable for space applications, but it makes it difficult to 
meet precise specified conditions. 

Project I was initiated not only to automate processing, but to guarantee that the 
specified mixing conditions will be precisely met. Temperature control of the 
propellants is critical and was taken as the first automation task. During mixing, a 
specified temperature profile must be met. 
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Project II will be based upon some of the existing knowledge and capabilities of 
computational simulation of engineering processes. The algorithms and simulation 
software that will be compiled during the course of this project will be provided to any 
research participant in the NASA Center. At the start of this project, the main objective 
has been focused on the computer simulation of the entire mixing process. 

Project III is concerned with the actual mechanical mixing of the propellants. 
Ultimately, one must be able to quantify the mixing process. Otherwise, how does one 
compare two mixtures? We must then seek out dynamical processes which yield 
"good" mixing in accordance with our definition. 

Approach and Results—Project I 

Two hot water tanks are currently in use with the University of Arizona mixer. One 
is set at 160°F and the other is set at 140°F. Currently, the operator draws water from 
one tank or the other to meet two temperature set points. This open-loop approach is 
undesirable, not only from the point of view that it requires an operator, but the 
temperature in the two tanks may vary considerably. Our approach is to continue to 
use the two tanks but with one set >160° and one set <140°. Two control valves will 
regulate the amount of water drawn from each tank, thereby maintaining a constant flow 
rate and allowing for control of the specified temperature set points. 

A proportional plus integral controller has been investigated using a numerical 
simulation and has been shown to be satisfactory. Actual control for the mixing system 
will be by means of a small digital computer with A/D and D/A interfaces. 
Thermocouple readings will be read into the computer and the control algorithm will, in 
turn, send a signal out through the D/A channel to set the valves. The computer system 
and one thermocouple are already in place. The control valves are on order and some 
replumbing will be required. 

Approach and Results—Project II 

Until we acquire a workstation for computer graphics and simulations, we are using 
the Silicon Graphics Workstation at the Computer-Aided Engineering Laboratory (CAEL). 
Since this is a Unix-based system, transition to another Unix-based system, such as a 
SUN workstation, should be rather easy. The first problem that we have considered for 
computer simulation is the formation of the mixing process of Project III. 

Approach and Results—Project III 

There is a close relationship between chaos and mixing. Indeed, the "baker map" 
used to explain chaos is a mixing tool (stretch then fold). This presents us with an 
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inverse problem in controls. Rather than designing a controller to produce uniform 
stable motion, our objective is to introduce controls into a system to make the behavior 
as random as possible. How does one design a controller so that the controlled system 
map will be one which will guarantee mixing? As with most progress in chaos theory, 
a great deal of simulation is required. We are currently working with some maps 
which should produce good mixing. These will be tested against a mixing definition so 
that a comparison of maps can be made. We will then examine the control dynamics 
to produce the maps. These will be examined both numerically and with an 
experimental mixing chamber. The first chamber has already been built and will be 
tested shortly. 
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Abstract 

Energy management opportunities in the process of hydrogen reduction of ilmenite 
for lunar oxygen production are being investigated. An optimal energy system to 
supply the power requirements for the process will be determined. 


Introduction 

Oxygen is one of the most important materials necessary for space exploration and 
colonization. It provides life support on space stations and fuel for propulsion. The 
high cost of payload delivery from Earth and the high fuel requirement make any long- 
distance space journey unfeasible. Oxygen production in space, such as on the moon 
or asteroids, can provide "refueling stations" for life support, combustion, and 
propellant. 

Different processes have been proposed for oxygen production in space. These 
processes will require high thermal and electrical power. This project examines the 
power requirements of one such process. Specifically, lunar oxygen production from 
ilmenite reaction with hydrogen will be studied to optimize its energy requirements. 


Objective 

The objective of the study is to determine an optimal energy system which would 
provide power for lunar oxygen production through hydrogen reduction of ilmenite. 
Ultimately, the goal is to develop an analytical model which would determine an 
appropriate optimal energy system for any given oxygen production process. 

Background 

The basic reactions for hydrogen reduction of ilmenite to produce oxygen are: 


FeTi0 3 + H 2 - *~ i 

Fe + Ti0 2 + H 2 0 

(1) 

h 2 o - 

H 2 + (1/2)0 2 

(2) 


The first reaction occurs endothermically at temperatures ranging from 700°C to 
1000°C. The per-pass conversion of H 2 to H 2 0 is about 5% at 900° C. The second 
reaction occurs through electrolysis, or thermal water splitting. As can be seen, there 
is a tremendous energy requirement, not just to meet the thermal and electrical 
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demands for the reactions to take place, but also for mining the ores, separating 
ilmenite from the ores, and pumping different materials through the system. 

Figure 3.2 illustrates the process with a simplified block diagram. Mined regolith 
first passes through a beneficiator, where ilmenite is separated from non-ferrous ores. 
The ilmenite is then sent to the reactor, where hydrogen is introduced. The hydrogen 
removes oxygen from the FeO component in ilmenite to form water. Water is then 
separated into components of hydrogen and oxygen through electrolysis. The oxygen 
is then cooled, liquefied, and stored. The hydrogen is recycled through the system. 

Approach 

The energy management study for the ilmenite hydrogen reduction process will 
include the following: 

1. Preparation of the ilmenite hydrogen reduction process flow diagram and 
determination of the appropriate temperatures and pressures for ilmenite feed, 
ilmenite/hydrogen reaction, water splitting, and the liquefied oxygen product. 
(This information will be developed with support from the Department of 
Chemical Engineering.) 

2. Use of computer programs (presently available or newly developed) to perform 
mass, energy, and exergy balance of all the components in the ilmenite 
hydrogen reduction process. This analysis will show the irreversibilities of the 
various components, identify the ones which are high contributors to 
irreversibilities, and, as a result, determine energy management opportunities. 

3. Application of the "Pinch Method" and use of commerically available computer 
programs to identify the various hot and cold streams in the process which 
could be used for process heat integration. Hot streams are those which need 
to be cooled; cold streams are those which need to be heated. The general 
approach of the Pinch Method is to combine all the hot streams in the process 
and plot them on an entha I py-versus- temperature chart. All the cold streams 
are combined and put on the same chart. The point where the smallest AT 
occurs beteen the two curves is the pinch. The hot streams above the pinch 
should be utilized for preheating purposes. The cold streams below the pinch 
should be utilized for precooling. Thus, the Pinch Method shows how various 
flow streams may be integrated to reduce heating and cooling requirements. 
For example, since it is desirable to reject heat in the spent feeds at as close 
to ambient temperature as possible, the high-temperature slag exiting the 
reactor may be used to preheat recycled hydrogen. 
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Figure 3.2. Block diagram for the process of hydrogen reduction 
produce oxygen. 


ilmenite to 









III-11 


4. Exploration of opportunities to minimize the use of electrical energy. For 
example, investigate the current status of research in thermal water splitting to 
determine if the technology could be applied to replace the electrolysis 
portion of the ilmenite hydrogen reduction process. 

5. Determination of the optimum source(s) of energy by setting the energy and 
capital (mass) targets. 

Research Status and Results to Date 

Due to problems with University registration, this project was not formally begun 
until January 1989 and, therefore, no analysis results have been obtained at this time. 

The research effort thus far has been limited to reading references on related 
subjects: 

1. W. W. Mendel I (Ed.), Lunar Bases and Space Activities of the 21st Century, 
Lunar and Planetary Institute, Houston, 1985. 

2. T. J. Kotas, The Exergy Method of Thermal Plant Analysis, Anchor Brendon 
Ltd., Essex, 1985. 

3. T. Rosenqvist, Principles of Extractive Metallurgy, McGraw-Hill Ind., 1974. 

4. Papers on the application of the Pinch Method to various industrial processes. 

Future Work 

Approach items (1) and (2) are expected to be completed during the next report 
period. The process flow diagram, as well as the process mass, energy, and exergy 
balance, will be presented in the next status report. The pinch analysis should also 
begin during this time. 
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Development of the data base for near-Earth resources was begun in the period 
covered by this report. 1 Specific items accomplished were: 

1. Definition of data base contents via discussions among the PSI staff and 
outside scientists such as Lucy McFadden and Mark Skyes. 

2. Discussion/demonstration of dBase capabilities and possible organization of 
file(s) with L. Alvarez, Steward Observatory, 

3. Data base entries for near-Earth asteroids consisting of data from the Russian 
ephemeris, a forthcoming review chapter for the Asteroids II book, and other 
literature sources. 

4. Development of a list of asteroid observers worldwide to be contacted during 
the coordination phase of the project. 

A copy of the data base as it now exists is included as Appendix C. It is 
currently set up in a spreadsheet format (which can easily be converted into an ASCII 
file), awaiting delivery of the dBase IV software from the University of Arizona as 
discussed with Dr. Lewis. We will transfer the data to dBase IV as soon as it arrives. 
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The Lunar and Planetary Institute (LPI) search produced approximately 80 
references. Most of these are available at the LPI. The product includes only -the 
author, title, date, and reference. The earliest (Cooper et al, Certain ecological aspects 
of a closed lunar base, Rand Corp.) is dated 6 March 1958. The latest is a semi- 
popular article (Cole & Majdacic, Astronomy, vol. 16, 18-24) produced last year. The 
early references tend to be reports of government-sponsored engineering studies. The 
later are generally semi-popular and science press articles. These references will be 
obtained on disk from the LPI and can be converted into an ASCII file for access via 
IBM-compatible PCs. 

At the Johnson Space Center (JSC) library a search was conducted of the major 
NASA data bases. The final print-out of references and short abstracts is 70 pages 
long. The document could not be obtained in electronic form through JSC. Most, 
likely all, of the useful references can be obtained through The University of Arizona 
library. Two days of travel and time to work at Tucson with the appropriate personnel 
will be allocated to obtain the desired references in computer-readable form. 

Twenty-three sets of microfiche were received from NASA on the earliest 
references to use of lunar materials. 

The Lunar and Planetary Receiving Lab (LPRL) at JSC has a small library with a 
very limited collection of materials on lunar utilization. However, an early study by Mr. 
A. Sanders (Extraterrestrial Consumables Production and Utilization, JSC, May 1972) was 
obtained and will be included. This report must be copied and the original returned to 
the JSC/Building 31 library. The University of Arizona library probably contains all the 
LPRL references of use to SERC. However, the group has produced "A Lunar 

Bibliography" (p. 50, March 1988) that should be of use. 

Following completion of the code Z Lunar Enterprise Study, the above materials 
and the reports on lunar processing done under my direction at the LPI will be 
reviewed for references of greatest interest to SERC. Senior researchers who 
participated in the workshops supported by that study will also be contacted to request 
copies of any relevant works they have produced since 1980. 
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An Annotated Bibliography of Propellant Processing Methods 

A. H. Cutler 

Dr. Cutler has specialized in space resources research for many years, with special 
emphasis on oxygen reduction methods. He has been retained to write a 
comprehensive review of these methods, detailing advantages and disadvantages, listing 
byproducts, and presenting recommendations. As the first step, an extensive outline 
was prepared, and the portion of this outline covering product manufacture was 
selected as the initial goal. A working version of this section fottows: 


Products 

I. Propellant 
A. Oxidizer 
1. Oxygen 

a. Production from lunar materials 

i. Slagging carbothermal 

- of ilmenite 

- of olivine 

- in integrated polyelement processes 

ii. Subsolidus reduction 

- ilmenite with H 2 or C 

- olinite with H 2 or C 

iii. Plasma processes 

- with H*/H 2 + FeO*X0 2 

- with CI-/CI 2 + FeO-X0 2 

iv. Slog electrolysis 

- bulk soil 

- mineral separates 

- with fluxes (Kesterk, EMEC) 

b. Production from asteroidal/Phobos/Deimos meet? 

i. Water splitting 

- electrolysis 

- photolysis 

- thermolysis 

ii. C0 2 processing 

- solid state electrolysis 

- carbonate? melt electrolysis 

iii. Oxide reductions, e.g., FeO + (CH 2 ) n 

c. Production on Mars surface— variants of processes above 
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2. Peroxide 

a. H 2 0 2 via Na 2 0 2 

b. Direct H 2 + 0 2 

c. Electrochemical means from 0 2 

d. Water partial oxidation 

3. Nitrogeneous oxidizers 

a. Organic N -* NH 3 

b. Organic N -*• N 2 

c. N 2 - NH 3 

d. Organic N -*■ NO 

e. NH 3 - NO + H 2 0 

f. NO - N0 2 

g. N0 2 + H 2 0 - HN0 3 

h. N 2 0 4 

i. RFNA 


B. Fuels 


a. Hydrogen 

i. From lunar soil 

ii. From asteroidal materials via hydrocarbons 

iii. Various routes from water (see above) 

b. Hydrocarbons 

i. Direct from asteroids phobos and demo 

ii. Reformed from initial pyrolyzate 

iii. Fischer-Tropsch synthesis from CO + H 2 0 

iv. From CH 3 OH 

c. Oxygenates— alcohols, ethers, etc. 

i. Methanol from CO + H 2 

ii. Higher alcohols from CO + H 2 

iii. Hydrocarbons ■+ oxygenates 

d. Hydrocarbon— hydrogen mixtures 
(solubility of H 2 in various materials) 

e. Metals 

i. Al 1 


Mg 


see metal production below 


iii. Ca J 

f. Metal/liquid mixtures 
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g. Nitrogeneous fuels 

i. Ammonia (see above) 

ii. Hydrazines from ammonia 

- Cl 2 process 

- 0 2 process 

- organo substituted hydrazines 

C. Thermal propellants 

1. Hydrogen (see above) 

2. Ammonia (see above) 

3. Hydrozine (see above) 

4. Others 

D. Elecric propellants 

1. Calcium and potassium 

2. Compounds— Cl 2 , S0 3 , etc. 

E. Mars propellants 

1. Hydrogen peroxide (see above) 

2. Hydrozine (see above) 

3. Fuel/oxidizer solutions 

a. Hydrocarbon— oxygen 

b. Peroxide/alcohol or peroxide/hydrocarbon 

4. Metal/sulfur propellants (MgS, CaS) 

II. Structural Materials 

A. Metallic 

1. Ferrous (Fe, Ni, Co) 

a. Fe— steel 

i. From magnetic collection 

ii. From subsolidus reduction 

iii. From slogging 

b. Iron alloys 

i. Imported alloying additions 

ii. Local source 

c. Forming and heat treating 

i. Forging 

ii. Casting 

iii. Machining 

iv. Carboxyl deposition 
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2. Non-ferrous metals 

a. Aluminum 

i. Winning 

- carbochlorination 

- plasma chlorination 

- fluorination 

- integrated processes 

ii. Forming 

- cast 

- roll 

- forge 

- machine 

iii. Alloying additions 

- local source— Mg, Mn 

- imported— Zn, Cu 

b. Magnesium 

i. Winning 

- carbothermal 

- salt electrolysis 

- plasma processing 

ii. Generation of MgO feed 

iii. Forming (as for aluminum) 

iv. Alloying additions 

3. Vacuum refining 

a. Advantages and effects on properties 

b. Effect on local environment 

4. Tool materials 

a. Tool steels— composition and properties 

b. Titanium carbide tool materials— cobalt bonded 

i. Cobalt source 

ii. Other binders 

5. Brazes and solders 

a. Aluminum, magnesium, and calcium as vacuum brazes 
B. Non-metal lie 
1. Bricks 

a. Morter bound 

b. Cut stone 
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c. Sintered 

d. Slogged and cast 

i. Glass 

ii. Ceramic 

2. Bogs 

a. Wire 

b. Glass 

3. Beams and pillars 

a. Fused soil 

b. Sintered 

c. Mortered 

d. Packed form 

4. Fibers— glass 

5. Refractories 

a. Graphite 

b. Magnetite and CaO 

c. CaS and MgS 

d. Cr 2 0 2 

e. Titanium oxides 

6. Cements 

7. Inorganic Polymers 

8. Sulfur 
III. Explosives 

A. Fuel/oxidizer explosives 

1. Fuel/H 2 0 2 

a. Organic 

b. Metallic 

2. NH4NO3 + fuel 

B. Nitrates from oxygenated organics and nitric acid 
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Preliminary Analysis of Surface Characteristics of Comet Tempel 2 

H. Campins 

Planetary Science Institute 

Science Applications International Corporation /,{ / 


A preliminary analysis of the surface characteristics of Comet Tempel 2 has been 
completed and incorporated into the paper "The Nucleus of Comet Tempel 2*" by M. F. 
A'Hearn, H. Campins, D. G. Schleicher, and R. L. Millis which was submitted to the 
Astro physical Journal, March 1989 (Appendix D). 

A detailed analysis of the surface characteristics of Comet Arend-Rigaux and 
Comet Tempel 2 is being carried out in collaboration with Dr. R. H. Brown (JPL). 


''Although the period of performance for Year 1 is given as 9/1/88-12/31/89, please 
note that funding (and therefore work) did not begin on these projects until 12/17/88. 
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Notes on Lunar Ilmenite 
M. Hutson 

Lunar and Planetary Laboratory 

The University of Arizona \ ^ 

Opaques (mostly ilmenite) make up 0-5% of highland rocks, 1-11% of low-Ti mare 
basalts, and 10-34% of high-Ti mare basalts (Carter 1988). Apollos 11 and 17 sampled 
high-Ti basalts. Apollos 12 and 14 sampled low-Ti basalts. Apollo 15 sampled a 
complex mixture of mare and highland material. Apollo 16 sampled mainly highland 
material (Taylor 1975). 



Mare Basalts (Rocks) 

Ilmenite is the third most abundant mineral in Apollos 11 and 17 rocks. It is very 
conspicuous in Apollo 11 samples, with platy crystals occasionally reaching several 
millimeters in size (Levinson and Taylor 1971, Mason and Melson 1970). Ilmenite is 
found in a wide variety of shapes, ranging from euhedral to anhedral, blocky, tabular, 
platy, and skeletal forms (Frondel 1975, Levinson 1971). Ilmenite also occurs as rims 
surrounding armalcolite cores. These rimming ilmenites generally contain exsolved 
rutile and spinel (Frondel 1975). However, most ilmenite is homogeneous, and the 
ilmenites associated with armalcolite and exsolved rutile are "relatively rare and 
quantitatively minor" (Levinson and Taylor 1971). 

Lunar ilmenite is almost pure FeTi0 3 ; a Mossbauer study of Apollo 11 ilmenite 
found no detectable ferric iron (Mason and Melson 1970). The following minor and 
trace elements have also been reported: Apollo 11— Al 2 0 3 , 0.1-0.3%; MnO, 0.3-0. 6%; 

CaO, up to 0.3%; Cr 2 0 3 , 0.1-1 .3%; V 2 0 3 , up to 0.2%; Zr 2 0 3 , up to 0.3%; and MgO, up 
to 6.0% (Mason and Melson 1970). Ilmenite appears to be the major host phase for Zr 
(Levinson and Taylor 1971). There does not appear to be any relationship between 
MgO content of ilmenite and its grain size or mineral association (Frondel 1975), 
although some researchers have suggested that high-magnesian ilmenite is associated 
with armalcolite (Frondel 1975). Levinson and Taylor (1971) state that high-magnesian 
ilmenite in Apollo 11 samples occurs as small tabular grains totally enclosed in 
pyroxene. Apollo 14 ilmenite contains 0.2-0.6 wt% Zr0 2 and up to 3.9 wt% MgO, which 
is noteworthy as the whole rock has lower MgO than that from other Apollo sites 
(approximately 8-8.5 wt%). In addition, ilmenite from one Apollo 14 sample contained 
Ti +3 (El Goresy et al. 1972). Apollo 17 sampled two types of basalts with different 
crystallization histories. Ilmenite is found as primary blocky crystals in both types of 
basalt and as rims around armalcolite in Type II basalts (El Goresy et al. 1974). Apollo 
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17 ilmenite contains 0.55-0.70 wt% MnO, 1.51-4.68 wt% MgO, and 0.48-1.05 wt% Cr 203 
(El Goresy et al. 1974). 

Mare Soils 

The lunar regolith contains material ranging from fines to blocks meters across. 
The material with grain size less than 1 mm is referred to as soil (Taylor 1975). This 
soil consists of mineral grains, rock fragments, glasses, and agglutinates, which are 
"glassy, rock and mineral fragments welded together by glass" (Taylor 1975). 
Agglutinates constitute up to half of the material in some size fractions (Levinson and 
Taylor 1971; Agosto 1988). As a result, the amount of ilmenite which is easily 
separated from lunar soils is considerably less than the normative ilmenite. This was 
shown by McKay and Williams (1979), who report modal abundances of 2.2% and 4.9% 
for Apollo 11 and 17 soils with normative abundances (calculated from Ti0 2 contents) 
of 14.7% and 15.5%, respectively. They note that the available ilmenite in the average 
Apollo 17 mare soil is only one-third of the normative ilmenite. McKay and Williams 
(1979) also suggest that ilmenite abundance increases with decreasing grain size 
fraction. However, a study of Luna 24 soils does not support this correlation (Laul et 
al. 1987). 

Frondel (1975) reports the results of x-ray studies; lunar and terrestrial ilmenite 
appear to have identical crystal structure parameters. She concludes: "It would appear 
that until we have many more x-ray and chemical analyses of both lunar and terrestrial 
ilmenites, the lunar ilmenite and its earthly counterpart can be assumed to be the 
same." However, it is known that reduction of lunar ilmenite is more difficult than for 
terrestrial ilmenite, due to the lack of ferric iron in the lunar ilmenite (Briggs 1988). 

Finally, McKay and Williams (1979) note that hydrogen is implanted in lunar soils 
by the solar wind. Data on hydrogen content of lunar soils are sparse and not in good 
agreeement. Reported bulk soil hydrogen contents range from 30-145 ppm. McKay 
and Williams (1979) note that ilmenite is known to retain helium much more readily than 
other minerals. Although there are no data for hydrogen, they speculate that if 
hydrogen is retained similarly to helium, it could exist at the 1000 ppm level in lunar 
ilmenite. 

Beneficiation 

A discussion of possible techniques for mining and beneficiation of lunar ores is 
given by Williams et al. (1979). However, it appears that only William N. Agosto has 
actually attempted to separate ilmenite from lunar soils in the laboratory (Agosto 1984, 
1985, 1988). Agosto has done electostatic separation experiments on simulant soils 
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(made with terrestrial ilmenites) and a 90-150 micron size fraction of an Apollo 11 soil. 
The experiments were run both in a vacuum and in a nitrogen atmosphere. Only 30% 
of the ilmenite was recovered in one pass in the vacuum (it is not clear whether he is 
referring to lunar or simulant ilmenites), while much better results were obtained using 
a nitrogen atmosphere, because of the density difference between ilmenite and other 
minerals in the soil (Agosto 1988). Ilmenite in the lunar simulant was concentrated from 
10 wt% to 95 wt%, with 68% of the ilmenite recovered after one pass in nitrogen. 
Ilmenite in the Apollo 11 soil was concentrated from 10 vol% to 54 vol%, with 45% 
recovery after one pass (Agosto 1984). Perhaps some combination of magnetic and 
electrostatic separation is needed. 
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ABSTRACT 


We present the results of simultaneous optical photometry and infrared ra- 
diometry of comet P/Tempel 2. Periodic variations of brightness, previously de- 
tected in the optical by others, are present at all wavelengths in our data. Because 
the optical and thermal lightcurves are in phase, we conclude that the variations 
are caused by the changing apparent cross-section of an elongated nucleus ro- 
tating with a period near 8.9 hours. The variation of flux with aperture allows 
us to separate the contributions of the nucleus and the coma. In the apertures 
with which we monitored the lightcurve, the contribution by the coma is roughly 
25% at maximum light in the optical and undetectable at the level of 10% at all 
times in the thermal infrared. By applying a standard thermal model, we have 
determined the following nuclear properties: effective radius at maximum light, 
5.9 km; visual geometric albedo, 0.022; and projected axial ratio near 2. The 
nucleus is very red. Based on the observed rate of outgassing, roughly 1% of the 
surface is active. Properties common to all nuclei of periodic comets appear to 
include the very low albedo and the high axial ratio. The active fraction of the 
surface and the spectral reflectivity appear to vary significandy among nuclei. 
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L INTRODUCTION 

Over the last ten years or so, a resurgence of interest in measuring the 
properties of cometary nuclei has occurred This resurgence is motivated by 
a desire to better understand the relationship between the nucleus and the more 
readily observed phenomena in the coma and to explore the relationship between 
cometary nuclei and Amor-Apollo (AA) asteroids. The comae of comets have 
been conceptualized for many years as the product of outgassing from a dirty 
snowball (Whipple, 1950) but the detailed nature of that snowball has remained 
elusive. Even the dramatic results from the Vega and Giotto spacecraft have left 
many questions unanswered concerning the nature of cometaiy nuclei, not the 
least of which is the degree to which P/Halley is typical of all comets. The 
possible relationship of comets to asteroids is based on the long-standing idea, 
originally suggested by Opik (1963), that a significant fraction of the AA asteroids 
are extinct cometary nuclei. Current estimates from dynamical arguments suggest, 
with very large error bars, that roughly half the AA asteroids are derived from 
comets (Wetherill, 1988). 

Unfortunately, it is very difficult to observe cometary nuclei. The only 
totally unambiguous observations are those from spacecraft, but ground-based 
observations have been made which, with varying degrees of certainty, appear to 
refer to a cometary nucleus rather than to material of the coma. These results 
have been reviewed recently (A’Heam, 1988; Weissman et aL, 1989) and will 
not be discussed in detail here. In general there are two possible approaches to 
the problem of studying cometary nuclei from Earth. Jewitt and Meech (1988) 
have concentrated on observing comets optically at large heliocentric distances 
when the activity of the comet is expected to be small. We, on the other hand, 
have concentrated on observing comets which are relatively inactive even at 
small heliocentric distances. This approach permits one to study the nuclei in 
the thermal infrared as well as in the optical (Campins et al., 1987; hereafter 
CAM87; Millis et al., 1983; hereafter MAC88). With both approaches there is 
always some uncertainty about the contamination of the measurements by coma 
although observers universally argue that their own measurements are relatively 
free from this problem. 

Comet P/Tempel 2 is a particularly interesting object. It has been considered 
repeatedly as a potential target for a cometary mission, but none of these has yet 
come to fruition. Nevertheless, because of this possibility, the comet has been 


l 
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the object of ground-based study by numerous astronomers. It is known from 
the orbital lightcurve (Sekanina, 1979) 1 that for much of the preperihelion arc 
the brightness of the coma is negligible. Similarly, the non-gravitational force is 
very small (Marsden, 1985) suggesting low activity per unit mass of the nucleus 
(although a fortuitous geometry could also lead to a small non-gravitational effect 
even with significant outgassing). P/Tempel 2 is, therefore, one of the few comets, 
the nuclei of which might be observable near perihelion. Accordingly, we planned 
an observational program of simultaneous optical and infrared photometry rimilar 
to those we carried out earlier on comets P/Neujmin 1 (CAM87) and P/Arend- 
Rigaux (MAC88). Other groups also observed P/Tempel 2 at this apparition. 
Jewitt and Luu (1988a) determined a rotational period from ccd photometry 
obtained well before perihelion, while Wisniewski (1988) confirmed the period 
and determined the rotational lightcurve from conventional broadband photometry 
in May 1988. Wisniewski’s rapidly distributed lightcurve was invaluable in the 
planning of our observing strategy. In particular, after the first few observations 
we were able to predict the times of subsequent maxima and minima in order to 
optimize the observations. 


H OBSERVATIONS 

The observations were scheduled for June 1988 near the time of the comet’s 
closest approach to Earth. Although this time was somewhat after the predicted 
onset of activity, the lightcurve of Sekanina (1979) implied that activity would 
still be weak and the advantage to be gained by observing at minimum geocentric 
distance was significant. The infrared observations were carried out at the NASA 
Infrared Telescope Facility ( IRTF) to measure the thermal emission from the 
comet on 1988 June 10 and 11 UT. The optical measurements were made with 
the University of Hawaii 88-inch telescope during the same two nights and during 
a small portion of the night of June 9 UT. Subsequent nights were cloudy. The 
observing circumstances at the midpoint of the two full nights are given in Table 1. 


1 An updated version of the orbital lightcurve, including data from the apparition of 1983, was 
widely distributed by the International Halley Watch and is the version actually used. 
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Tabic 1 

Observational Circumstances June 10.9 


rn = 1.72 AU phase = 20° 

A = 0.79 AU = —9.9 km/s 


In the infrared, the standard N filter (reference wavelength = 10.1 /xm; 
Tokunaga, 1986) and IRTF bolometer were used most frequently to define the 
lightcurve. The temporal variations were monitored using the 4— mm aperture. 
A few N observations were made near maximum and minimum light (phases 
chosen to minimize the temporal variations) in 4 apertures, 3, 4, 5, and 6.5 
mm, for which drift scans on stars yielded effective angular diameters of 5.3, 
7.0, 7.8 and 8.7 arcsec, respectively. The chopper throw was 30 arcsec north- 
south. Based on the photometry with different apertures, we estimate that the 
coma contribution to the reference beam was less than 0.5% of the comet’s peak 
brightness. Measurements with the standard M (4.8 p m) and Q (20.0 /im) filters 
were also made to determine color temperatures. The star a Bootis was used as 
a primary standard; the adopted magnitudes were M = -2.12, N = -3.17 and 
Q = —3.13 (Tokunaga 1986). Two stars near the comet in the sky, a Serpentis 
and 110 Virginis, were observed frequently and used as local standards to monitor 
extinction and check for any instrumental variations. Neither of the local standards 
showed intrinsic variability although some changes in extinction were detected and 
corrected for. The observations were reduced and calibrated using the procedure 
and absolute fluxes discussed by Tokunaga (1986). 

The optical observations were made with a photometer and pulse-counting 
electronics from Lowell Observatory and the standard International Halley Watch 
(IHW) filters which isolate various emission bands and portions of the reflected 
continuum. These filters are described in detail in an appendix to the IHW 
Archive of Observations of PlGiacobini-Zinner (in press). On June 9 an EMI 
6256S phototube was used completing the identical system to that which was used 
previously for comet P/Arend-Rigaux (MAC88). On June 10 and 11, we used an 
RCA 31034A phototube which allowed use of filters at longer wavelengths than 
were used in the previous study. Since most of June 9 was devoted to another 
project, only about 1 1/2 hours were spent observing PyTempel 2. Consequently, 
the data from that night were not well calibrated. The temporal variations were 
monitored using the two filters for the continuum at 4845 A and 6840 A and 
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an aperture of 10.1 arc sec. Sky was observed far from the comet (usually 
several tens of arcmin) after each set of observations of the comet (usually 
3 thirty-second integrations in each of the two continuum filters). The stars 
HD 149363 and HD 120086 were used for absolute calibration via the standard 
magnitudes defined by the IHW while observations of two nearby comparison 
stars were interspersed between every set of cometary observations to monitor 
any variations in transparency and/or instrumental sensitivity. There was no 
indication of intrinsic variability in the brightness of the nearby comparison stars 
or in the extinction coefficients. Mean extinction coefficients for the two nights 
were used in reducing the data. Two or three times each night — usually near 
the times of maximum and minimum light when the brightness of the comet was 
changing most slowly — observations were made in apertures of 14.2, 20.0, and 
28.5 arcsec in order to determine the radial brightness profile of the coma. An 
accurate measurement of this profile is required if the contributions of the nucleus 
and the coma to the observed signal are to be disentangled. The comet was also 
observed a few times through other filters of the IHW set in order to estimate 
production rates of various gasses, gas-to-dust ratio, etc. 

A. Brightness Variations 

The 6840-A observations taken through the 10.1-arcsec aperture and the N 
magnitudes measured through a 7— arcsec aperture are listed in Tables 2 and 3 and 
are plotted as a function of time in Figures 1 and 2 for the nights of 10 and 1 1 
June, respectively. It is evident that both the reflected optical and the thermal 
infrared brightnesses of the comet were varying in a seemingly periodic way. 
The existence of this variation was, in fact, no surprise. Jewitt and Luu (1988a) 
had reported cyclic variations in April with an approximately 9-hour period and a 
peak-to-peak amplitude of 0.7 mag. in the R passband, in good agreement with the 
amplitude we observed at N, while Wisniewski (1988), from V-filter observations 
in late May, had produced a convincing light curve for Tempel 2 having a period 
of 8 hr. 58 min. and an amplitude of 0.5 mag. Wisniewski’s observations were 
made through photometer apertures of 17 and 12 arc sec diameter; Jewitt and 
Luu observed with a CCD. The difference in optical amplitude reported by these 
investigators and the yet smaller amplitude of approximately 0.4 mag which we 
observed in June at 6840A, very probably can be attributed to greater dilution 
by the coma of the signal from the nucleus as the comet became more active 
because, based on Sekanina’s (1987) determination of the polar orientation, we 
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all observed the comet at nearly identical projections. On the other hand, the 
observations of P/Tempel 2 in 1987 by Jewitt and Meech (1988), which exhibited 
a much smaller amplitude (<0.3 mag given our present knowledge of the period), 
were obtained when the sub-Earth and sub-solar points were at comctographic 
latitudes corresponding to a projection factor approximately half that for the later 
observations. Consequently, the smaller optical amplitude in 1987 is likely due 
to the different aspect. 







■ | 


Although our observations do not extend over a sufficient time interval to 
permit a precise period determination, we have estimated the period using the 
method of phase dispersion minimization (Stcllingwcrf 1978). This technique 
was previously applied to comet P/Arend-Rigaux by MAC88 and the reader is 
referred to that paper for details of the method. We performed separate period 
searches on the 4845-A, 6840-A, and 10-pim data for P/Tempel 2 from 10 and 
11 June. The statistic 9 (see MAC88) for each of these data sets is shown plotted 
against period in Figure 3. Possible periods in the data are those for which there 
is a local minimum in this plot The increase in depth of the minima towards 
shorter periods is, to some degree, an artifact of the method. In computing 9, 
the data are divided into phase bins and the dispersion of the data points within 
each bin computed. The shorter the period, the narrower is each bin in time, 
and the less the intrinsic brightness variation of the comet will contribute to the 
dispersion within a bin. Note that the plots for all 3 data sets show a sharp 
minimum at a period of 4.45 days. This period corresponds to a light curve 
having one maximum and one minimum per cycle. We believe, however, based 
on evidence which will be presented later in this paper, that the light curve must 
have 2 maxima and 2 minima per cycle. Consequently, the true period is twice 
the value quoted above, in good agreement with the results of Jewitt and Luu 
(1988a, 1989) and Wisniewski (1988). Indeed, a broad minimum is present in 
Figure 3 at a period near 9 hours. We note that this is the shortest photometrically 
determined rotational period for a cometary nucleus except for the very first period 
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so determined. Fay and Wisniewski (1978) determined a period of 5.2 hours for 
P/d’ Arrest although the rather small amplitude of the variation yielded poor signal- 
to-noise ratio in the resultant power spectrum. All other rotational periods recently 
determined seem to be substantially longer. 

Figure 4 shows the observations from 10 and 11 June phased according to a 
period of 8.9 hours and plotted as open circles and filled circles, respectively. In 
addition, the observations from 9 June have been included (as filled triangles) in 
the plot of the 4845-A data. Zero phase arbitrarily has been set at 7.75 hrs. U.T. 
on 9 June. Since the observations from 9 June are poorly calibrated, we have, 
simply shifted those data arbitrarily to the same mean brightness as the data at the 
same phase from the other nights. Note that the apparent peak in the abbreviated 
light curve from 9 June at phase 0.1 is coincident with a maximum in the light 
curve from the other two nights. 

It is apparent in Figures 1, 2, and 4 that the reflected optical and thermal 
infrared light curves are closely in phase as was also the case for P/Arend— 
Rigaux (MAC88). Now, as in the earlier paper, we interpret the alignment of 
the visible and infrared light curves to mean that the brightness variations result 
from the changing apparent cross section of an elongated, rotating nucleus. If the 
brightness variations were due to albedo features on the nucleus, the visible and 
IR light curves should be 180° out of phase. The large difference in amplitude 
between the visible and infrared light curves results at least in part from greater 
dilution of the optical observations by the signal from the coma. In order to 
derive the size, shape, and albedo of P/Tempel 2’s nucleus, it is necessary to first 
remove this contamination by the coma. 

B. Removal of Coma 

As noted above, observations were made of the aperture-dependence of the 
flux several times at the three monitoring wavelengths — 4845 A, 6840 A, and 
10.1 fim. These observations were usually made near maximum or minimum 
light in order to minimize the temporal variation of the comet’s brightness 
during the interval required to make measurements through all apertures. The 
sequence of these measurements was such that an observation with the monitoring 
aperture (10.1 arcsec and 7.0 arc sec in the optical and infrared, respectively) was 
interpolated between each of the observations with the other apertures. This 
procedure took up to an hour so it was necessary to use smooth curves through 
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the data in Rgures 1 and 2 to correct the observations to a common dme assuming 
that the variation in flux was the same in all apertures, an assumption justified a 
posteriori. These corrections were as large as 6%. 

A simple radial outflow model for the dust predicts a surface brightness which 
varies as p~ ! , where p is the projected distance from the nucleus, which in tum 
implies that the flux in a centered aperture varies linearly with the diameter. 
Although this is clearly an oversimplification since the coma may be varying, it 
provides a convenient reference point. In Rgures 5 and 6 we have plotted the 
flux as a function of diameter of the aperture. An examination of the variation 
with aperture size at 10.1 shows that a horizontal straight line adequately 
describes the data. Rtting straight lines by least squares and comparing the value 
of those lines at 7.0 arcsec and 0 arcsec shows that the contribution by the coma 
to the monitoring aperture is 10% + 10% at maximum light Thus, there is not 
a statistically significant detection of the coma at this wavelength and we assume 
that the contribution is nil for our reductions. There was not enough signal to even 
attempt to measure the variation with aperture size with the M and Q filters and 
therefore we have assumed that the contribution by the coma at those wavelengths 
is also negligible. At least for Q, this assumption is supported a posteriori by 
the modelling discussed below. 

At optical wavelengths, however, there is clearly a contribution by the coma. 
In Figure 6, the straight lines have been fitted by least squares to observations 
through the 10.1, 14.2, 20.0, and 28.5-arcsec apertures. The standard deviations of 
the slopes are roughly 10% and the variation in slope from one set of observations 
to the next is not statistically significant. Based on the average slopes from Figure 
6, the contribution by the coma to the flux in the 10.1-arcsec aperture at maximum 
light is 25% of the total flux at 4845 A and 20% of the total flux at 6840 A. At 
minimum light, the contributions are 35% and 28%, respectively. 

The lack of variation in the optical brightness profiles, it should be noted, 
argues strongly that the coma is not contributing significantly to the observed 
brightness variations. For one thing, the amplitude of the observed variations is 
markedly greater than the variation in the coma contribution to the signal in the 
10.1-arcsec aperture. Secondly, since the crossing time for dust in our apertures is 
comparable to the rotational period, any substantial variation in coma production 
would produce readily apparent distortions of the coma brightness profiles. 

The coma contribution to the observations of the comet in the 10.1-arcsec 
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aperture has been calculated on the basis of these values of the slope and subtracted 
from the observed signal to give the brightness of the nucleus alone. Columns 3 
and 7 of Table 2 contain these nuclear magnitudes on the standard system defined 
by the IHW . These values arc plotted as a function of lightcurve phase in Figure 
7. The observed magnitudes at N, given in Table 3, are also plotted in Figure 
7. We have already argued that these 10.1— /rm observations directly refer to the 
nucleus and do not require correction for coma. 

A careful comparison of the panels in Figure 7 shows that the amplitude of 
the optical lightcurves is about 0.55 mag. while that in the infrared is near 0.80 
mag. There are two possible explanations of this difference. One is that we have 
underestimated the amount of coma contamination of the optical measurements. In 
fact examination of Figure 6 does suggest a tendency for the observed variation 
with aperture to deviate from linearity in the sense of being convex upwards. 
Curvature of this type has been seen in some other comets and has been attributed 
to the effects of radiation pressure (Jewitt and Meech, 1987) or evaporation of 
grains (Baum and Kreidl, 1986). To explore the impact of possible curvature in 
the brightness profile of the coma, quadratic fits to the data in Figure 6 were 
performed. These fits yielded nuclear fluxes at maximum light 10 to 50% smaller 
than deduced from the linear fits and, therefore, gave lightcurve amplitudes closer 
to that observed in the infrared. However, a reduced chi— squared analysis indicates 
that the quality of the quadratic fits is not significantly improved over that of 
the linear fits. Accordingly, we see no compelling justification for adopting the 
quadratic fits in computing the contribution by the coma. This is particularly true 
in the absence of any physical understanding of why the brightness profile should 
follow a quadratic function. 

An alternative explanation of the difference in amplitude between the optical 
and IR lightcurves comes from the work of Brown (1985). Brown has shown 
that for elongated bodies, the amplitude of the thermal signal will be a function 
of wavelength as well as aspect. Because we have only the 10.1— ^m lightcurve 
for P /Tempel 2, we are unable to confirm whether this effect is actually present in 
this comet. If it is, then it is by no means clear that the optical and IR lightcurves 
should have the same amplitude. In any case, it is important to emphasize that if 
the coma contribution has been underestimated, the result will be an overestimate 
of the albedo of the nucleus but, because the albedo turns out to be so low anyway, 
the determination of the nuclear size will be unaff ected 
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C. Emission by the Coma 

The observations of the coma, aimed at deriving rates of outgassing of the var- 
ious molecular species, were also interspersed with the monitoring observations. 
These observations were, therefore, much less complete than if determination of 
the outgassing were our primary goal The only species measured with large fields 
of view were OH- and CN, the former because it is the best measure of the total 
gaseous outflow and the latter because it is usually the easiest species to measure 
in comets. C2 and C3 were measured only in the 10.1 arcsec aperture used for 
monitoring the variability. Of these, C2 was measured numerous times because, 
if strong, it could contaminate the measurements of the continuum at 4845 A. 
Fortunately, the emission was weak enough that the contamination of the 4845-A 
flux in the 10.1-arcsec aperture was negligible although still strong enough that 
we were able to detect it reliably with the filter centered on the emission band 
(5140 A). 

The observations of OH required special treatment because of the nonlinearity 
of atmospheric attenuation with airmass in the filter used to isolate this band. The 
IHW has modelled the atmospheric attenuation in this filter (Carsenty, private 
communication) using the actual transmission curve of this filter and atmospheric 
parameters suitable for Mauna Kea. The extinction can be linearized for these 
specific conditions by using an effective airmass: 

X eff = X- 0.07257 x X 2 ; (B - stars ) 

X e ff = X — 0.0348 x X 2 \ ( comets ) 
k comet — 0-9234 X kg 

Conversion to fluxes, both for OH and for the other species, was carried out using 
the transformations provided by the IHW ( cf. appendices to the Giacobini-Zinner 
archive, in press). These formulae require continuum measurements which are 
used to remove the continuum that underlies the emission bands. In all cases the 
measurements of the continuum were sufficiently close in time that no correction 
was needed for the rotational variation. The resultant fluxes in the emission 
bands are given in Table 5 below. 

The dust in the coma was measured as a byproduct of the removal of the 
brightness of the coma from the nuclear brightness as described above. This 
process is subject to all the uncertainties discussed in the previous section to an 
even higher degree since the desired component is now the smaller of the two 
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components being separated- A single measurement at 4845 A in an aperture of 
56.8 arcsec is less susceptible to these problems (only 28% of the light was due 
to the nucleus) but even in this case the result is somewhat uncertain. 


HI. MODELLING AND INTERPRETATION 


A. Coma 


We interpret the observed emission-band fluxes in the coma in terms of 
molecular production rates as is common in cometary studies. The fluorescence 
efficiencies (g-factors) used to convert fluxes to column densities are given in 
Table 4. The column densities are then converted to production rates using a 
Haser model with the parameters given in Table 4. From the abun dan ces of OH, 
we also calculate the production rate of H 2 O using a vectorial model (Festou, 
1981) with the parameters given in Table 4. For the species for which we have 
multiple measurements, there is no statistically significant variation with time so 
we have averaged all measurements to give an average production rate for June 
10-11 for a given aperture. These results are given in Table 5. An examination 
of the relative abundances of the various gaseous species shows that the ratio of 
other species to H 2 O is somewhat lower than typical of active comets but not 
by a large factor. 
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Table 4 

Haser Model Parameters 
(at r=1.72AU) 


Species 

UN 

[erg-s -1 ] 

Parent Scale 
[km] 

Daughter Scale 
[km] 

OH 

8-23X10 -16 

1.20x10 s 

3.41x10 s 

CN 

1.16xl0' 13 

3.60x10* 

8.83 x10 s 

c 2 

1.53 xlO- 13 

4.71x10* 

3.24x10 s 

c 3 

3.39xl0' 13 

2.94 x10 s 

1.77 x10 s 


Vectorial Model Parameters 



(at r=l AU) 


r(H 2 0) 

= 82400 s 

r(OH) 

= 200,000 s 

v(OH) = 

1.15 km/s 

fraction H 2 0 to OH = 0.93 


The production of dust is given by the quantity Af p (A ’Hearn et al ., 1983), 
a directly observable (i.e., model independent) quantity which, for pure radial 
outflow at constant velocity and for constant scattering, is proportional to the 
production of those grains dominating the optical scattering. Since it was shown 
above that the coma was constant to within our precision, we quote only a single 
value for each aperture. A few measurements were made with a filter at 3650 A 
(see Table 7 below) through the 10.1-arcsec aperture but, as with M and Q, there 
was not sufficient signal-to-noise ratio to separate the coma from the nucleus 
directly. On the other hand, these data clearly show that the amplitude of the 
lightcurve at this wavelength is roughly 0.7 magnitudes, noticeably greater than 
that found at longer optical wavelengths. These data will be discussed in more 
detail below but, since all evidence points to a constant coma, the amplitude 
suggests that the contribution by the coma is less than 20%, thereby allowing us 
to set an upper limit on the brightness of the coma at this wavelength. 

We note that all our modelling is based on models that assume spherical 
symmetry whereas it is well know that the coma of P/Tempel 2 is not circularly 
symmetric ( e.g . Sekanina, 1987). In many cases the lack of spherical symmetry 
will not materially affect the conclusions since, e.g., a set of jets in the plane of 
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the sky can be modelled as simple fractions of the spherically symmetric case. 
Deviations from spherical symmetry become systematically less important as the 
field of view increases and in the limiting case in which one observes the entire 
comet the discrepancies vanish. For this reason we place much heavier weight on 
the observations with the large (56. 8-arcsec) aperture. Temporal variability is not 
included in any of our models and one must remember that any ’production rates’ 
represent an average rate over some appropriate interval depending on the transit 
time of material in the field of view and the lifetimes of the relevant species and 
their parents. 


Table 5 

Production of the Coma 
Gas 


Species 

Aperture 

Radius 

Flux 

N 

Q 


Diam. 


[erg-cm' 2 - 




[arcsec] 

[10 4 km] 

s' 1 ] 

[cm" 2 ] 

[s- 1 ] 

OH 

56.8 

1.62 

3.3 xlO' 12 

8.4xlO n 

2.1 xlO 27 


10.1 

0.29 

1.5 xlO" 13 

1.17xl0 12 

1.9 xlO 27 

H 2 0 

56.8 

1.62 

- 

- 

3.2X10 27 


10.1 

0.29 

- 

- 

3.8 xlO 27 

CN 

56.8 

1.62 

1.6xl0* 12 

4xl0 9 

3x 10 24 

c 2 

10.1 

0.29 

5 x 10" 14 

8 x10 s 

2x 10 24 


10.1 

0.29 

3xl0‘ 14 

8 x10 s 

6x 10 22 



Dust 



Wavelength 

Aperture 

Radius 

Flux 

Alp 


Diam. 





[A] 

[arcsec ] 

[10 4 km] 

[erg-cm" 2 - 

■s-'-A- 1 ] 

[cm] 

3650 

10.1 

0.29 

<3xl0' 16 

<18 

4845 

56.8 

1.62 

3.0xl0" 15 

15 


10.1 

0.29 

6.1 xlO" 16 

18 

6840 

10.1 

0.29 

4.4 x 10" 16 

17 
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The ratio of gas to dust is conveniently parameterized by the ratio Q(H20)/Afp. 
This ratio is independent of the effects due to the finite field of view of the instru- 
ments although it differs from the true ratio of gas to dust both by a conversion of 
units (Afp has the units of cm) and by parameters of the dust and the gas which 
may vary from time to time and from comet to comet In particular, it varies 
with the wavelength of the light used to determine Afp due to the variation of 
reflectivity with wavelength. Nevertheless, we believe that this is the most obser- 
vationally reasonable parameter for monitoring the ratio of gas to dust Feldman 
and A’Heam (1985) showed that the parameter did exhibit a correlation with 
heliocentric distance, which might indicate a variation of modelling parameters 
but might also indicate real physical differences with heliocentric distance. At a 
given distance, the comets considered there showed rather little scatter. Comet 
P/Halley exhibited significant variation with heliocentric distance and closer than 
1 AU did deviate noticeably from the usual relationship with distance (Feldman 
era /., 1987). Typical values of this ratio in P/Halley are in the range 10 25 to 10 26 
(cm-s) -1 using ultraviolet light and these numbers would be very roughly 2/3 that 
value if the continuum were measured at 4845A. Comets generally considered to 
be dust-poor like P/Encke and Bradfield 1979 X have values near 2x 10 26 using 
the continuum at 4845 A. The value for P/Tempel 2 is roughly lx 10 26 [cm-s] -1 , 
a very typical value for comets at this heliocentric distance. The only comet of 
which we are aware that exhibits a substantially higher value of this parameter 
is P/Neujmin 1 (CAM87) for which the dust was undetectable and for which we 
estimate a lower limit of lxlO 27 [cm-s] -1 . 

B. Nucleus 

The nucleus was modeled using the "standard" thermal modeling program de- 
scribed by Lebofsky et al. (1986) as in the previous modeling of comets P/Neujmin 
1 (CAM 87) and P/Arend-Rigaux (MAC 88). The input parameters used were 
the following: thermal phase coefficient 0.01 mag/degree, visible phase coeffi- 
cient 0.035 mag/degree, thermal emissivity 0.9, and thermal beaming factor 0.9. 
The visible phase coefficient is consistent with but not strongly constrained by 
the data at phases < 30° presented by Sekanina (1987, Figure 11). None of the 
other coefficients are known for this or for any comet but are taken from the 
standard model for asteroids. Because the observations were made at a relatively 
small phase angle (Table 1), the results of the model are not very sensitive to as- 
sumptions about phase coefficients. This type of thermal modeling is appropriate 
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for comets that are relatively inactive, such as P/Tempel 2, because the amount 
of energy that goes into sublimation is negligible compared to the total incident 
radiation (CAM 87). The calculated effective radius at maximum light, which 
depends only on the infrared measurements as long as the optical albedo is low, 
is 5.90 +0.24/-0.68 km. In this calculation a 10.1-pm magnitude at maximum 
light of 3.532 was used. This value includes the monochromatic correction of 
-0.043 magnitudes for 270K, given by Tokunaga (1986). The larger uncertainty 
towards smaller radii is due to the possibility that the coma contributes, as much 
as 20% of the 10.1-pm flux at maximum light (see section 2). 

The axial ratios obtained from the corrected visible lightcurve amplitude (0.55 
mag) and the 10.1pm lightcurve amplitude (0.S0 mag) were averaged to determine 
a minimum axial ratio for the nucleus of 1.9 +0.6/-0.2; the larger uncertainty 
towards a more elongated nucleus is again due to the possibility of a larger 
contribution by the coma. At the time of our observations. Earth was nearly in 
the equatorial plane of the comet according to the model of Sekanina (1987); 
if so, and assuming that the nucleus of this comet can be approximated by a 
prolate spheroid rotating about a minor axis, the observed lightcurve amplitude 
and effective radius are consistent with overall nuclear dimensions of 16.3 x 8.6 
x 8.6 km. The geometric albedo is 0.022 with an estimated uncertainty of +.004/- 
.006. The effects of departures from sphericity on the results of the standard 
thermal model have been discussed by Brown (1985). We estimate the error in 
the calculated effective radius introduced by ignoring sphericity to be +/- 5%. 
Brown’s non-spherical thermal model predicts that the amplitude of the rotational 
lightcurve will be larger in the thermal infrared than in the visible. This effect 
could explain the different amplitudes observed in the visible and at 10.1 pm; 
on the other hand, this difference could also be due to an incomplete removal 
of the coma contribution in the visible {cf. section II.). A more refined analysis 
of our observations and those of MAC 88 aimed at improving our knowledge of 
the shape of the nuclei of comets P/Tempel 2 and P/Arend-Rigaux is in progress 
(Brown et al., in preparation). 

The observations at M (4.8 pm) and Q (20 pm), presented in Table 6, provide a 
check on the thermal modeling. Several measurements at Q were made at different 
points in the lightcurve in search of possible temperature variations predicted by 
the non-spherical thermal model. We detected only statistically marginal N— Q 
color variations, a mean N— Q = 1.85 was observed, consistent with the brightness 
temperatures of 270 K at N and 256 K at Q, predicted by the standard thermal 
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model. The comet was considerably fainter at M, hence the fluxes in this bandpass 
were not as well determined as those for Q. Only single values for M were 
obtained for each day near a maximum in the lighteurve; these are slightly higher 
than predicted by the standard thermal model. This excess in the M flux is 
qualitatively similar to that observed in comet P/Arend-Rigaux (MAC 88), where 
a significant contribution (approx. 33% of the flux) from the dust was observed 
in the M bandpass. In the case of comet PyTempel 2, however, we could not 
obtain aperture photometry in the M filter to quantify the contribution by the dust 
coma at this wavelength. 


Table 6 

Additional Thermal Data 


UT 

[hr] 

M 

[mag] 

UT 

[hr] 

Q 

[mag] 

June 10 


June 10 


12.48 

8.56 

13.05 

1.98 

June 11 


June 11 


10.25 

8.28 

07.60 

2.34 



08.98 

1.99 



09.63 

1.73 



10.72 

1.94 



11.58 

2.13 


The albedo of the nucleus derived here is an effective visual albedo but some 
additional information is available on the variation of albedo with wavelength. The 
narrowband magnitudes given in Table 2 yield directly the color of the nucleus 
(averaging values near maximum light), Cl (48 - 68) = m (4845) - m (6840) = 
1.18 whereas the solar color in this system is 0.87 based on observations of 
solar analogs. The nucleus of comet P/Tempel 2 is therefore rather red, with a 
reflectivity gradient, (ds/cLA)/s = 14% per kA. This result is in excellent agreement 
with the spectrophotometry of Tempel 2 at large heliocentric distance obtained 
by Spinrad et ai, (1979) and in reasonably good agreement with the results 
contemporaneous with ours by Jewitt and Luu (1989). The three sets of data arc 
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shown in Figure 8. Since the value of the albedo-area product by Spinrad et al. 
is also in good agreement with our value near minimum light, we infer that no 
coma was present at the time of the measurements by Spinrad et al.. Those data 
suggest that the reflectivity increases at the shortest wavelengths. 

Although we have data at 3650 A, we can not reliably remove the coma 
from those data because the signal-to-noise ratio is inadequate. Our data for 
the combined nucleus and coma at four different times are given in Table 7. 
Uncertainties arc roughly 0.02 mag at the two longer wavelengths and rather larger 
but not well determined at 3650 A. The table shows quite clearly the redness 
of P/Tempcl 2 between 4845 A and 6840 A with only marginally significant 
variations between maximum and minimum. The color index between 3650 A 
and 4845 A, however, varies from bluer than the sun at maximum light to slightly 
redder at minimum light. Since all data at longer wavelengths imply that the 
coma is constant, and since the amplitude at 3650 A is 0.7 mag (averaging the 
two maxima and comparing with the single minimum), significantly greater than 
at longer optical wavelengths, we conclude that the contribution by the coma 
at maximum light is less than 20%. We can therefore set a lower limit to 
the nuclear brightness: m nttc (3650) < 17.1 at maximum light yielding a color 
Cl (36 — 48) = m(3650) — m(4845) < 1.0 compared to the corresponding 
solar color of 1.17. This corresponds to a slope of the reflectivity of — 13 %/kA. 
The agreement among the three sets of data in Figure 8 is excellent and the 
behavior at short wavelengths, although based on only two points with limited 
precision, is remarkable. The two points at the shortest wavelengths in Figure 8, 
one of ours and one from Spinrad et al., are individually no more than three 
standard deviations above an extrapolation of the data at longer wavelengths but 
the agreement between the two totally independent results is striking. Given the 
uncertainties in the data, we do not want to overinterpret this result but believe 
that the effect is real. The repeatability of the spectrum over many different 
aspects and, at least for the filter photometry, with rotational phase indicates that 
the surface is rather homogeneous in color. 
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Tabic 7 

Color of P/Tempel 2 
(10.1 arc see aperture) 


Date UT 

m(3650) 

[mag] 

m(4845) 

[mag] 

0(36-48) 

[mag] 

0(48-68) 

[mag] 

Jun 10.475 

16.97 

15.90 

1.08 

1.20 

Jun 11.333 

17.56 

16.23 

1.34 

1.16 

Jun 11.404 

16.73 

15.82 

0.91 

1.08 

Jun 11.471 

17.36 

15.97 

1.40 

1.12 

solar color 

— 

— 

1.17 

0.87 


It is also significant that the color of the nucleus is quite different from the 
color of the grains in the coma. Examination of Table 5 shows that the grains in 
the coma have a very neutral color whereas the nucleus is distinctly red except at 
the shortest wavelengths. Jewitt and Luu (1989) have reached a similar conclusion 
based on their long-slit spectra of the comet There is no way of determining, 
from either our data or those of Jewitt and Luu, whether the difference is due to 
selective ejection of grains or to differences in scattering properties between the 
grains on the surface and grains in the optically thin coma. 

We can combine the outgassing rate observed in the coma and the size of the 
nucleus to estimate the fraction of active area on the nucleus. Determination of 
that fraction depends, of course, on where the active area is on the nucleus since 
that determines the effective insolation. We have used the methods described by 
Cowan and A’Heam (1979) to estimate the outgassing of water molecules per unit 
area for an isothermal nucleus (rapid rotadon and high thermal conductivity), for 
a model which is everywhere in local equilibrium (low thermal inertia) and for a 
model in which parallels of latitude are isotherms (rapid rotation). The range of 
rates of outgassing is from 2x 10 16 to 5x 10 17 cm -2 s -1 . According to the nuclear 
model of Sekanina (1987), the vent is at a high cometographic latitude. For this 
latitude and obliquity 90° and each latitude line being an isotherm, a typical rate 
of outgassing would be of order 3xl0 16 cm -2 s -1 . This is a lower limit to the 
outgassing per unit area since concentration of the heating to the daytime side 
non-linearly increases the outassing so that the average over the day is higher. 
Thus if the outgassing is from a single active vent vaporizing in equilibrium, its 


17 





D-28 


maximum area is 11 km 2 or 3% of the total surface. Much smaller areas are 
likely since the outgassing per unit area increases if the active area is more nearly 
normal to the direction to the sun or if a latitude line is much hotter on the dayside 
than on the nightside as we would expect It is clear from this discussion that 
comet P/Tempel 2, like comets P/Arend-Rigaux and P/Neujmin 1, has an almost 
totally inactive nucleus. The fraction of the surface which can be active in all 
three comets is at least an order of magnitude less than the fraction for P/Halley 
which itself has activity from no more than 30% of the surface (based on 15% in 
the sunlit hemisphere and an assumed comparable area on the dark hemisphere). 

We also note that the geometry of Sekanina’s model predicts that the vent 
should be alternating between daylight and darkness at the time of our observations 
so that some variability might be expected (Sekanina, 1988). Because the lifetimes 
of H 2 O and OH are both longer than the rotational period of the nucleus, the 
observed variation will be much less than the true variation in outgassing. A 
variation of the outgassing from one day to the next has been reported (A’Heam 
et al., 1988) and will be discussed in more detail elsewhere. On the other hand, 
the constancy of the coma contributed by the dust argues against significant 
fluctuations in the outgassing unless the velocity of the dust leads to a transit 
time in our aperture much longer than the rotational period. If the velocity of 
the dust was <0.1 km/s this would be the case and we should not see significant 
fluctuations. 


Table 8 

Summary of Nuclear Properties 

Reff = 5.9 km Pv = 0.022 Axial Ratio >1.9 

P = 8.9 hr Tb = 270 K Active fraction <3% 

prolate spheroid; 16.3x8.6x8.6 km 


The combination of the red slope of the reflectivity at long optical wavelengths 
with a possible increase in the ultraviolet is unusual. The only similar curve of 
spectral reflectivity for a small body of which we are aware is that of comet 
P/Schwas smann-Wachmann 1 (S-Wl). Although we can not be sure that the 
quiescent phase of this comet is a bare nucleus, that point is often advocated and 
not unreasonable. Cochran et al. (1982) measured the reflectivity of S-Wl in its 
quiescent phase and found a reflectivity similar to that of P/Tempel 2. Also as for 
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P/Tempel 2, they found that the coma of S-Wl (measured during a subsequent 
outburst) was of a different color. The two comets are not, however, identical 
since the albedo of S-Wl was found by Cruikshank and Brown (1983) to be 0.13, 
very much higher than found here for P/Tempel 2 

The reflectivities of these two comets appear unlike those of any standard 
types of asteroids. They are also significantly different from the reflectivities of 
other recently measured cometary nuclei: P/Neujmin 1 (CAM 87) was very red 
even at the shortest wavelengths; P/Arend-Rigaux (MAC88) was almost neutral 
in color, P/Halley (Thomas and Keller, 1989) was significandy less red at optical 
wavelengths. These other comets are similar in visual albedo and axial ratio and 
show a range of active surface fractions spanning the fraction found above for 
P/Tempel 2. The present results for P/Tempel 2 thus support the contention that 
cometary nuclei form a diverse group with no clear-cut match among the AA 
asteroids. 

We note finally that the results found here bear on models of cometary 
nuclei. The results discussed here are generally consistent with the model of 
Sekanina (1987) in regard to the polar orientation and the general geometric 
description of the outgassing as being from a vent at high cometographic latitude. 
The combination of these results with similar results on the fraction of the 
surface which is active on other cometary nuclei qualitatively supports the widely 
discussed picture involving a mantle which gradually builds up on the surface of 
the nucleus, thereby choking off the vaporization. 

IV. SUMMARY 

Simultaneous photometry at infrared and optical wavelengths has allowed 
us to determine the albedo, size, and projected shape of the nucleus of comet 
P/Tempel 2. In gross size and shape it appears to be nearly identical to the 
nucleus of P/Halley, a prolate spheroid with axes 16x8 1/2x8 1/2 km. It differs 
from that of P/Halley in that the active portion of the surface is much less, of order 
1% as opposed to an estimated 15-30% for P/Halley. The rotational period of 
the nucleus is near 9 hours, the second shortest of the photometrically determined 
rotational periods for comets. The outgassing and the changes in the observed 
amplitude of the rotational modulation appear to be consistent with the nuclear 
model for this comet by Sekanina (1987) although other aspects of the model (such 
as the assumption that the nucleus is oblate) appear to need revision. The spectral 
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reflectivity of the nucleus is red in the optical, but there is a suggestion that the 
reflectivity increases in the ultraviolet. In the optical the redness is significantly 
greater than for P/Halley and comparable to the redness of some but not all other 
cometary nuclei and D asteroids. The increase at ultraviolet wavelengths, if real, 
is similar only to that seen for P/S-Wl. 

Observations of the nuclei of several periodic comets in recent years have 
made it clear that they are all dark and usually very aspherical. Furthermore, 
most of the surface is inactive on the nuclei of all periodic comets for which 
the relevant data exist There is a selection effect in this conclusion because 
nuclei active over their entire surface are much harder to detect inside the 
coma. Nevertheless, a picture of gradual development of a mande culminating in 
extinction of activity qualitatively fits all the available data on cometary nuclei. 
It is still premature to make statements regarding the distribution of sizes or 
rotational periods because the sample of objects is too small and, at least as far 
as sizes are concerned, subject to very strong selection effects. Finally, it seems 
clear that nuclei differ significandy from one another in their spectral reflectivity. 
Although there arc many similarities to D and C asteroids, there also appear to be 
significant differences. For this reason, there is as yet no physical connection that 
can be direcdy established between cometary nuclei and Amor-Apollo asteroids. 
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FIGURE CAPTIONS 

Figure 1. Observed lightcurve of comet P/Tempel 2 in reflected optical light and 
thermal infrared radiation on 1988 June 10 UT. The optical observations were 
made with a 10.1— arcsec aperture and the infrared observations were made with 
a 7.0-arc sec aperture. 

Figure 2. Same as Figure 1 but for 1988 June 11. 

Figure 3. Phase -dispersion-measure for the lightcurve of comet P/Tempel 2 as a 
function of period. The statistic is shown for each of the three wavelengths at 
which lightcurves were obtained, 4845 A, 6840 A, and 10.1 /im. All show their 
deepest minimum at 4.45 hours which yields a single-peaked light curve and we 
believe that the correct rotational period is the one indicated by the broad but 
shallower minimum at twice that period, 8.9 hours. 

Figure 4. Phased lightcurves for comet P/Tempel 2 assuming a period of 8.9 hours 
including data from 1988 June 9 (4845 A only; filled triangles), June 10 (open 
circles), and June 11 (filled circles). Zero phase is at 1988 June 9, 7.75 hrs UT. 
Figure 5. Variation of flux with aperture at 10.1 ptm. The data are consistent 
with the absence of any coma. 

Figure 6. Variation of flux with aperture in the optical region. Straight lines have 
been fitted by least squares and indicate that a coma is definitely present. The 
average slope was used for removal of coma from the data. 

Figure 7. Rotational lightcurve of the nucleus of comet P/Tempel 2. The optical 
observations have been corrected for the contribution by the coma. The differ- 
ence in amplitude between the optical and infrared could be due to incomplete 
correction for coma but could also be an effect due to the enhanced amplitude 
of a prolate spheroid in its thermal emission. Zero phase is 1988 June 9, 7.75 
hours UT. 

Figure 8. Spectral reflectivity of the nucleus of comet P/Tempel 2. The two points 
at the shortest wavelengths differ from the extrapolated relationship from longer 
wavelengths by only about three standard deviations but appear to represent a 
real feature. 
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Abstract 

Carbonaceous chondrites are usually 
believed to be the primary constituents of 
near-Earth asteroids and Phobos and 
Diemos, and are potential resources of 
fuels which may be exploited for future 
planetary missions* In this work we have 
calculated the nature and abundances of 
the major volatile bearing and other 
phases, including the vapor phase, that 
should form in Cl and C2 type carbonaceous 
chondrites as functions of pressure and 
temperature. The results suggest that 
talc, antigorite ± magnesite are the major 
volatile bearing phases and are stable 
below 400°C at 1 bar in these chondritic 
compositions. Simulated heating of a 
kilogram of C2 chondrite at fixed bulk 
composition between 400 and 800°c at 1 bar 
yields about 135 gm of volatile, which is 
made primarily of H^O, CH 4 , C0 2 and 
CO. The relative abundances of these 
volatile species change as functions of 
temperature, and on a molar basis, H 2 
becomes the most dominant species above 
550°C. In contrast, Cl chondrites yield 
about 306 gm of volatile under the same 
condition, which consist almost completely 
of 60 wt% H 2 0 and 40 wt% C0 2 . Preliminary 
kinetic considerations suggest that equi- 
librium dehydration of hydrous phyllosili- 
cates should be attainable within a few 
hours at 600°c. These results provide the 
framework for further analyses of the 
volatile and economic resource potentials 
of carbonaceous chondrites. 

Introduction 

The near-Earth asteroids are poten- 
tial resources of volatiles which can be 
used as propellants and life supporting 
purposes for planetary missions. As 
emphasized by Lewis and Lewis (1), these 
asteroids are at times the nearest bodies 
to Earth, and many of them can pass 
between Earth and Moon. These authors 
have also shown that within any arbitrary 
three-year period, there are roughly 300 
and 90,000 launch opportunities to respec- 
tively kilometer and 100-meter sized 
asteroidal bodies as compared to 39 lunar 
launch windows. These facts and the 
potential volatile resources make the 
asteroids important candidates for 
detailed analyses for supporting future 
planetary missions. 

Cost-effective or energy-efficient 
extraction of volatiles from the asteroids 
require knowledge of the nature and modal 


abundances of the constituent minerals in 
which the volatile components are struc- 
turally bound. There has been several 
studies on reflectance spectroscopy of 
near-Earth asteroids to characterize the 
nature of the volatile bearing minerals 
(2,3). While these studies report absorp- 
tion band near 3 urn suggesting presence of 
hydrous phyllosilicates, more precise 
characterizations were equivocal. An 
obvious alternative approach to the reso- 
lution of the problem is to examine 
directly the asteroidal materials. 

The spectral and density character- 
istics strongly suggest that a significant 
fraction of the near-Earth asteroids are 
made of carbonaceous chondrites (1,4). 
The most volatile rich of these meteorites 
are what are commonly known as Cl or (Cl) 
and C2 (or CM) classes. However, the 
grain size of the volatile bearigg phases 
are often too small (100 to 1000 A) (5) to 
permit identification under optical micro- 
scope, and are best observed under high 
resolution transmission electron micro- 
scopy ( HRTEM ) . Zolensky and McSween (5) 
have recently presented a summary of vola- 
tile bearing phases identified in the 
carbonacous chondrites. 

The primary objective of this work is 
to carry out theoretical calculations to 
predict modal abundances and compositions 
of the major mineral phases, along with 
the abundance and composition of the co- 
existing vapor phase, that could develop 
in the bulk compositions of Cl and C2 
chondrites as functions of pressure and 
temperature. The results would provide 
the framework for engineering designs for 
the extraction of volatile components from 
asteroids as well as Phobos and Deimos, 
the two small natural satellites of Mars, 
which are also likely to be made of car- 
bonaceous chondrites (1). 

Theoretical Method 

Principles 

According to Duhem's theorem in 
classical thermodynamics (6,7), the equi- 
librium state of a closed system (i.e., a 
system of fixed composition and mass) is 
completely determined if any two variables 
are fixed regardless of whether these are 
intensive, extensive, or a combination of 
both. Ganguly and Saxena (7) have recent- 
ly reviewed the various methods by which 
one may carry out the actual computation 
of the equilibrium assemblages. The 



Table 1 Bulk chemical 
compositions of Cl and C2 
carbonaceous chondrites 
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Cl 

C2 


Si 

10.40 

12.96 

Ti 

0.04 

0.06 

Al 

0.84 

1.17 

Cr 

0.23 

0.29 

Fe 

18.67 

21.56 

Mn 

0.17 

0 . 16 

Mg 

9.60 

11.72 

Ca 

1.01 

1.32 

Na 

0.55 

0.42 

K 

0.05 

0.06 

P 

0.14 

0.13 

Ni 

1.03 

1.25 

Co 

0.05 

0.06 

S 

5.92 

3.38 

H 

2.08 

1.42 

C 

3.61 

2.30 

0 

45.61 

41.74 


method used in this study obtains the com- 
positions and abundances of the equilib- 
rium phase assemblages by minimizing the 
Gibbs Free Energy (G) of the system at 
fixed P-T conditions. The minimization is 
constrained to conserve the various ele- 
mental masses of the system through the 
method of Lagrangian multipliers (7). 
Because of the extremely small grain size 
of the phyllosil icates observed in carbo- 
naceous chondrites, one should consider 
the effects of surface free energy on the 
formation of these phases. As a first 
approximation, we have, however, ignored 
this effect. Therefore, the actual equi- 
librium temperatures of volatile bearing 
phases could be lower than those predicted 
by our calculations. 

System, Phases and Data Base 

Table 1 shows the bulk compositions 
of Cl and C2 carbonaceous chondrites. 
These compositions are based on analyses 
of a number of fragments of Orgueil and 


Table 2 Phases in the system 
Mg-Fe-Si-C-H-O-S considered in this work 

The sources of thermochemical data are shown by reference 
numbers within parentheses. The Heat of Formation of the Fe-end 
members of the hydrous phyllosilicates are derived in this work. 


Anhydrous Silicates and Oxides 


Ol ivine 

{ Mg , Fe ) 2 Si0 4 

( 8 ) 

Orthopyroxene 

(Mg,Fe)Si0 3 

( 8 ) 

Per iclase 

(Mg,Fe)0 

(23) 

Quartz 

Si0 2 

(17) 

Cristobal ite 

Si0 2 

(17) 

Hematite 

Fe 2°3 

(23) 

Magnetite 

Fe 3 °4 

(23) 

Hydrous and Carbonate 

Phases 


Anthophy 1 1 i te 

(Mg,Fe) 7 Si 8 0 2 2<OH) 2 

(8,27) 

Talc 

(Mg,Fe) 3 Si 4 O 10 <OH) 

( 8 ) 

Ant igori te 

( Mg , Fe ) 43^34035 (OH ) 82 

( 8 ) 

Chrysot ile 

( Mg , Fe ) 3 S i 2 0 $ ( OH ) 4 

( 8 ) 

Bruc i te 

(Mg , Fe ) ( OH ) 2 

(17,24) 

Magnesite 

{ Mg , Fe ) CO 3 

(17,25) 

Sulfides and Elements 



Troili te 

FeS 

(17) 

Pyr i te 

FeS 2 

(17) 

Iron 

Fe 

(26) 

Sulfur 

S 

(17) 

Graphite 

c 

( 17) 

Vapor Phase 

C-O-H-S 

(15,16,17) 

Species: H 2 Q, C0 2 , CH 4 , CO, 0 2 , H 2 S, S 2 , S0 2 , COS 
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Murchison meteorites, respectively, and 
are summarized in Dodd (22). It can be 
easily seen from this table that the sub- 
system Mg-Fe-Si-C-H-O-S (MFCHOS) consti- 
tutes almost 96% by weight of the bulk 
compositions of Cl and C2 chondrites. In 
order to somewhat simplify the computa- 
tional problem, we have confined the G 
minimization calculation to this sub- 
system, but evidently the phases whose 
compositions form outside of this sub- 
system cannot be major constituents of 
these chondrites. 

Table 2 shows the list of phases that 
we have considered as possible crystal- 
lizing phases within the above subsystem 
along with the sources of the thermochem- 
ical data. The choice of these phases 
have been guided by the published reports 
of the mineralogy of carbonaceous chon- 
drites whose compositions lie within the 
system MFCHOS (e.g., Zolensky and McSween 
(5)). The selection of thermochemical 
data in the system Mg-Si-O-H is based on 
the critical analysis of Chatterjee (8,11) 
for consistency with both calorimetric and 
experimentally determined phase equilib- 
rium measurements. These data are similar 
to those of Berman (9). However, unlike 
those in Berman, the data are also consis- 
tent with high P-T experimental data 
involving distribution of Fe and Mg among 
coexisting silicates. 

The heat capacities (C p ) and entro- 
pies (S°) of the phases, if *not available 
in any of the recent systematizations of 
self-consistent data set, are estimated by 
stoichiometric summation of the properties 
of oxides and structurally analogous com- 
pounds, as discussed in a number of publi- 
cations (13,14) • 

There are no thermochemical data for 
the Fe-end members of the phyllosilicates 
listed in Table 2. The enthalpies of 
these phases are estimated according to 
the known relative properties of Fe- and 
Mg-end members of biotite, which is also a 
phyllosilicate. The data derived by 
Chatterjee (8) and Robie 2 |nd Hemingway 
(10) for the Mg- and Fe^ -end members 
(phlogopite and annite, respectively) 
yield an enthalpy difference of -359.566 
KJ per mole of divalent cation at 1 bar, 
298 K. We have used these values to esti- 
mate the enthalpies of formation from 
elements of Fe-end members of phyllosili- 
cates at 1 bar, 298 K from the available 
data for the Mg-end members. The results 
are as follows (KJ/mol of Fe). Fe-talc: 
-1606.80; Fe-antigorite: -1127.53 and Fe- 
chrysotile: -1094.87. The corresponding 
values of the free energy of formation 
from elements (A^f) for Fe-talc and Fe- 
antigorite at 1 bar, 29^ K are -4477.18 
and -4784.0 KJ/mol of Fe^ , respectively. 

The validity. of the estimated thermo- 
chemical properties of the Fe-end members 
of the phyllosilicates may be tested as 
follows. The Fe-end members of talc and 


antigorite are unstable. Instead, the 
stable phases are closely analogous com- 
pounds, namely, minnesotite, which has a 
stoichiometry of Fe 27 S i 36°86 ( 0H ) 26 as 
compared to that of Fe27Si 38 09Q ( OH ) jo for 
Fe-talc, and greenalite with a stoicniom- 
etry of Fe4 8 Si 32 0 8 Q < 0H > 64 compared to 
Fe 48 Si 3 40 8 5(0H) 62 for Fe-antigorite. 
Anovitz et al. (12) have recently deter- 
mined AGf of minnesotite and greenalite 
from phase equilibrium data. Their 
results, AG° (minnesotite) = -4577.78 KJ 
and AGf (greenalite) = -4799.50 KJ, 
differ in the right direction from the 
corresponding values of Fe-talc and Fe- 
antigorite, considering the relative 
stabilities of the phases. 

The stereochemical environments of 
divalent cations in biotite are quite 
similar to those in talc, both being 2:1 
(TOT) layer silicates, and partly similar 
to those in antigorite, which is an 1:1 
(TO) layer silicate (28). Thus, these 
phyllosilicates should have similar 
differences between the enthalpies of 
formation of the Fe- and Mg-end members. 
Further, since on the basis of observa- 
tional data, the predicted compositions of 
the mineral phases in carbonaceous chon- 
drites are expected to be Mg-rich, the 
effects of errors in the Fe-end member 
thermochemical properties and solid solu- 
tion model will be relatively small. This 
can be easily understood by considering 
the form of G-X curve of a stable solution 
and evaluating the effects of errors in 
above properties near Mg-terminal segment. 

The ±hermodynamic mixing properties 
of the Fe^ - and Mg-end member components 
of the anhydrous silicates are taken from 
the self-consistent summary of Chatterjee 
(8), while those of the phyllosilicates 
are assumed to be ideal for the lack of 
any data. The latter assumption is likgly 
to be approximately valid at T _>_ 600 C, 
since at these conditions biotite solid 
solution is found to behave approximately 
ideally (7). At lower temperature, the 
phyllosilicates may deviate from ideal 
solution behavior, but perhaps not too 
dissimilarly to significantly affect their 
relative stabilities. Magnesite is also 
assumed to be an ideal solution of Fe-Mg 
components. As emphasized by Ganguly and 
Saxena (7), the computed dehydration con- 
ditions are not very sensitive to errors 
in activity-composition relations since 
the devolatilization equilibria are char- 
acterized by large enthalpy changes. The 
iron oxides and sulfides are considered to 
be stoichiometric phases since the solid 
solutions of Mg are around a few percent 
at the temperatures of interest. 

The fluid properties are taken from 
Saxena and Fei (15,16) and Robie et al. 
(17). The "correspond ing state method" of 
estimating P-V-T relation, as discussed 
by these authors, is extended to include 
the sulfur species. For mixing data, we 
have assumed that H 2 S behave as H 2 0, ^2 as 
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Fig, 1 Calculated modal abundances of minerals per kilogram of 
rock in (a) Cl and (b) C2 carbonaceous chondritic bulk compositions 
as functions of temperature at 2 Kb, The steps in the calculations 
are shown by symbols. 


C> 2 , and SO 2 and COS as CO 2 # and used the 
equivalent binary interaction parameters 
from Saxena and Fei (15), 

Results 

At 1 bar, the free energy minimiza- 
tion calculations do not yield any vola- 
tile bearing phase at T >_ 400 C. At the 
present state of the thermochemical data, 
we do not believe that our calculations 


are reliable at lower temperatures. In- 
deed, we find results at lower tempera- 
tures which do not seem to follow the 
systematics of the results obtained at 
higher temperatures. Thus, we have calcu- 
lated the equilibrium phase assemblages at 
a pressure of 2 Kb so that the volatile 
bearing phases appear at higher tempera- 
tures. This procedure at least permits 
recognition of the major volatile bearing 
phases. Some of these phases may- be re- 
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Fig. 2 Calculated compositions of the silicate phases in Cl and 
C2 carbonaceous chondrites as functions of temperature at 2 Kb, 


placed by more stable assemblages at lower 
temperature, but should reappear during 
heating experiment. The results, which 
are illustrated in Figure 1, show that the 
primary volatile bearing phases in the Cl 
and C2 bulk compositions are talc, anti- 
gorite ± magnesite. Talc is stable to 
around 650°C at 2 Kb, and dehydrates to 
orthopyroxene (OPx) at higher temperature 
according to Tc + 01 * 5 OPx + H 2 O. Anti- 
gorite and magnesite crystallize at T < 
550°c, but magnesite is absent in C2 bulk 
composition at T 2. 400°C. FeS is a stable 
phase at all temperatures investigated in 
this work. We, however, did not consider 
the formation of sulfate or pyrrhotite. 
Magnetite becomes a stable phase at 
T 500°C. Elemental sulfur found in Cl 
chondrite is invariably associated with 
pyrrhotite (Fei_ x S) (5). This suggests 
decomposition of troilite (FeS) to pyrrho- 
tite and sulfur. 

Figure 2 shows the compositions of 
the f erromagnes i an phases as function of 
temperature at 2 Kb for the Cl and C2 bulk 
compositions. Calculations at 1 bar also 
yield similar compositions. The observed 
compositions of these phases in the car- 
bonaceous chondrites vary widely in Fe/Mg 
ratio. One way to achieve this diverse 
Fe/Mg ratio is by aqueous alteration which 
seems to be a characteristic of the car- 
bonaceous chondrites. Fe z fractionates 
very strongly into an aqueous phase rela- 
tive to the silicates (18). Thus, varying 
degrees of equilibration with an aqueous 
phase could lead to a spectrum of Fe/Mg 
compositions of the ferromagnesian sili- 


cates. Aqueous alteration may not be the 
sole reason, but is likely to be an impor- 
tant contributing factor in the develop- 
ment of variable Fe/Mg ratio of the ferro- 
magnesian silicates in carbonaceous chon- 
drites. 

The abundances of volatile species in 
equilibrium with the solids per kilogram 
of total mass of C2 chondrites are illus- 
trated in Figure 3. Figures 3a and 3b 
illustrate the molar and mass abundances, 
respectively, as functions of temperature 
at 1 bar. It should be noted that the 
total mass of the volatiles do not change 
above 400°C, which means that under equi- 
librium condition C2 chondritic material 
will completely devolatize if heated above 
400°c. The total yield of the volatile is 
about 14% of the initial mass of C2 chon- 
drites. The change in th^ abundance of 
various species above 400 C is due to 
homogeneous reactions within the fluid 
phase. Above 550 C, hydrogen has the 
highest molar concentration in the vapor 
phase. Hydrogen is a very good propellant 
and reducing agent which may thus be pref- 
erentially extracted by heating the C2 
chondritic material to T > 550 c. 

In contrast to the C2 chondrites, the 
volatiles given off by Cl material above 
400°C consist essentially of H 2 0 and CO 2 
at 1 bar pressure. Per kilogram of mate- 
rial, the total yield of volatile is about 
306 grams of which nearly 60 wt% is H 2 0. 

The above results are subject to 
considerable uncertainties owing to the 



E-6 




Fig* 3 Calculated equilibrium (a) molar and (b) mass abundances of 
the volatile species per kilogram of C2 chondritic material. 


approximate nature of the thermochemical 
properties of the phyllosil icate solid 
solutions, and also since we have ignored 
the effects of surface free energy on the 
formation of extremely fine grained crys- 
tals* Further, carbonaceous chondrites of 
any given type do not have homogeneous 
compositions, but instead probably repre- 
sent agglomeration of materials formed at 
different conditions at different parts of 
the solar system. However, these calcula- 
tions are still useful in providing an 
idea of the relative stabilities, as well 


as upper stabilities of the volatile bear- 
ing phases and their average abundance. 
It is interesting to note that recent 
reflectance spectroscopic studies of car- 
bonaceous chondrites seem to corroborate 
the predictions made above as to the 
nature of the most abundant hydrous phases 
(Gaffey, pers* comm.). 

Kinetics of Devolatilization 

The equilibrium calculations pre- 
sented above provide a framework for the 
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7Talc - 3Anth + 4Qz + 4H20 



Fig. 4 Dehydration of talc as functions of time, temper- 
ature and extent of overstepping above the equilibrium 
dehydration boundary at 1 Kb. Volume fraction is that of 
talc remaining after reaction. Experimental data points, 
shown by asterisks, are from Greenwood (19). 
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Fig. 5 Equilibrium relations of talc, 
antigorite and brucite along with other 
minerals in the system MgO-S iC^-^O. 
Reproduced from Evans and Guggenheim (21). 


evaluation of the volatile as well as some 
minera log ical resource potential of Cl and 
C2 carbonaceous chondrites. The actual 
extraction of volatiles as a function of 
temperature and heating rate, however, 
depends on the kinetics of the devolatili- 
zation process. For kinetic analysis, it 
is essential to identify the minerals in 
which the volatiles are most likely to be 
structurally bound. The above results 
suggest that talc, antigorite and magne- 
site are likely to be the major volatile 
bearing phases in carbonaceous chondrites. 

Except for pure Mg-talc, there are no 
kinetic data for the other volatile bear- 
ing phases. The dehydration kinetics of 
talc was determined by Greenwood (19). At 
a given temperature, the rate o| reaction 
depends on the rate constant, K r and the 
extent of departure from the equilibrium 
boundary. This is expressed in transition 
state theory (20) according to the rela- 
tion Rate * K ( 1-exp ( nAG/RT ) ) , where n is 
a constant, commonly assumed to be unity, 
and AG is the Gibbs free energy of the 
reaction at the P,T condition of interest. 
Figure 4, which is constructed from Green- 
wood's experimental data, illustrates the 
dehydration kinetics of talc as functions 
of time, temperature and extent of over- 
stepping (AT) above the equilibrium dehy- 
dration boundary at 1 Kb. It is evident 
from this figure that essentially complete 
dehydration of talc can be achieved, within 
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6 hours at about 160°C above the equilib- 
rium boundary. 


The stability fields of the pure Mg- 
end members of various phyllosi 1 ica tes of 
interest in this study are summarized for 
the condition of Ptptal 3 P H20 F ^9 ure 
5, which is reproduced from Evans and 
Guggenheim (21). According to our re- 
sults, talc and other phyllos il ica tes com- 
pletely dehydrate at >_ 400°C at 1 bar in 
Cl and C2 bulk compositions. This is com- 
patible with the phase relations shown in 
Figure 5. Thus, based on the results 
shown in Figure 4, one would anticipate 
that essentially complete dehydration of 
all phyllosi 1 icates in Cl and C2 chondrit- 
ic material is likely to be achieved ig a 
matter of a few hours at T _>_ 600 c. 
Detailed studies of the devolatilization 
kinetics of both phyl losi 1 icates and car- 
bonates as functions of temperature, fluid 
composition, and grain size are currently 
under way to evaluate the practical limi- 
tations of the extraction of volatiles 
from carbonaceous chondritic materials. 
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